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I.  INTRODUCTION 


Knowledge  of  the  makeup  of  the  Earth’s  atmosphere  is  an  area  of  intense  academic 
and  government  interest.  Specifically,  the  Department  of  Defense  (DoD)  has  a  stated 
requirement  to  gain  an  understanding  of  the  electron  density  profile  of  the  ionosphere  for 
applications  in  over-the-horizon  RADAR  and  long-range  communications.  Knowledge  of 
the  altitude  distribution  of  atomic  oxygen  is  of  interest  to  the  scientific  conununity 
because  atomic  oxygen  density  is  known  to  be  an  impoitant  component  of  photochemical 
models  of  the  atmosphere.  Additionally,  DoD  and  NASA  have  become  concerned  with  the 
material  degradation  of  spacecraft  surfaces  exposed  to  atomic  oxygen  in  flight.  In  each 
application,  accurate  measurement  of  the  altitude  distribution  of  atomic  oxygen  is  a 
requirement. 

The  most  common  method  for  studying  the  upper  atmosphere  has  been  through  the 
use  of  experimental  sounding  rockets.  Additional  information  has  been  gained  through 
extensive  study  of  the  propagation  of  electromagnetic  waves  in  the  ionosphere.  The  Naval 
Postgraduate  School  Physics  Department,  through  the  efforts  of  the  Center  for  Space  Sys¬ 
tems  Applications,  has  been  involved  in  ionospheric  research  for  several  years.  They  have 
conducted  two  successful  sounding  rocket  experiments  in  conjunction  with  the  Naval 
Research  Laboratory  (NRL)  and  are  funded  for  a  third  experiment  in  eariy  1994.  The  NPS 
portion  of  these  experiments  has  been  passive  sensing  of  the  middle  ultraviolet  region  of 
the  electromagnetic  spectrum  by  the  Middle  Ultraviolet  Spectrograph  (MUSTANG) 
instrument. 

The  evolution  of  this  research  has  brought  forth  a  desire  to  examine  the  Far  and 
Extreme  Ultraviolet  regions  of  the  electromagnetic  spectrum  (FUV  and  EUV,  respec¬ 
tively)  at  moderate  to  high  resolution.  Moderate  resolution  instruments  allow  the  study  of 
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atmosphere  dynamics  through  the  examination  of  the  Doppler  shifts  of  emission  lines. 
High  resolution  instruments  have  the  ci^>ability  to  measure  emission  line  profiles  and  thus 
pennit  the  investigation  of  radiation  transport  through  the  examination  of  optically  thick 
lines.  (Cleary,  et  al.,  1992) 

A.  THESIS  OBJECTIVES 

The  main  objective  of  this  thesis  project  was  to  validate  the  design  of  a  high  resolu¬ 
tion  All-reflection  Michelson  Interferometer  for  use  in  the  FUV  region  of  the  electromag¬ 
netic  spectrum.  This  validation  was  to  be  accomplished  through  the  use  of  computer  ray 
tracing  procedures  and  fourier  analysis.  A  secondary  objective  was  to  further  examine  and 
optimize  various  parameters  of  the  instrument  in  order  to  propose  design  revisions.  The 
interferometer  design  was  originally  examined  in  detail  by  Wallace  (1992)  and  Professor 
David  D.  Cleary  of  the  NPS  Physics  department.  This  research  is  a  continuation  of  those 
efforts. 

B.  THESIS  OUTLINE 

Chq>ter  n  of  this  thesis  provides  background  information  regarding  the  composition 
of  the  atmosphere.  The  importance  of  atmospheric  research  and  the  significance  of  this 
project  in  the  area  of  high  resolution  spectroscopy  are  also  discussed.  The  third  chr^ter 
gives  a  basic  review  of  the  theory  of  interference  and  di^action  and  an  introduction  to  the 
All-reflection  Michelson  Interferometer  (AMI).  The  simulation  of  the  Oxygen  1304  A 
triplet  emission  line  profile  is  then  presented.  Chapter  IV  introduces  the  ray  tracing  soft¬ 
ware  used  and  discusses  its  limitations.  The  fifth  cht^ter  explains  the  procedure  used  in 
conducting  this  research  and  presents  the  specific  results  obtained.  Finally,  Ch{q>ter  VI 
offers  conclusions  and  makes  recommendations  for  further  study  in  this  area.  Two  qTpen- 
dices  are  included  which  contain  ray  tracing  parameter  files  and  program  listings. 
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II.  BACKGROUND 


A.  ATMOSPHERIC  COMPOSITION 

The  Earth’s  atmosphere  is  generally  divided  into  regions  based  on  kinetic  tempera¬ 
ture  or  on  chemical  con^sition.  The  four  main  temperature  regions  or  “sj^res”  of  the 
earth  have  characteristic  temperature  gradients  and  are  separated  by  nearly  isothermal 
transition  regions  referred  to  as  “pauses”.  The  lowest  temperature  shell  is  the  troposphere, 
the  area  of  the  atmosphere  where  surface  weather  exists.  The  troposphere  has  a  negative 
ten^rature  gradient  and  extends  from  the  Earth’s  surface  to  an  altitude  of  approximately 
ten  kilometers.  It  is  followed  at  an  altitude  of  10-12  km  by  the  tropopause.  The  second 
region,  extending  from  approximately  12  km  to  SO  km  in  altitude,  is  the  stratosphere.  The 
stratosphere  maintains  a  positive  temperature  gradient  mainly  due  to  heating  caused  by  the 
ozone  layer  (discussed  below).  It  is  followed  by  the  stratopause.  Above  the  stratopause  is 
the  mesosphere,  characterized  by  a  negative  temperature  gradient.  Above  the  mesosphere 
is  the  mesopause,  followed  by  the  thermosphere  and  thermopause.  The  thermosphere  is 
heated  primarily  by  the  absorption  of  electromagnetic  radiation  by  molecular  oxygen  in 
the  range  of  1000  A  to  2000  A.  This  wanning  initiates  a  final  temperature  gradient  rever¬ 
sal.  In  the  thermosphere,  interaction  between  solar  radiation  and  atmospheric  gases  forms 
the  region  of  ionized  particles  known  uS  the  ionosphere.  The  temperature  regimes  and  the 
associated  kinetic  temperature  profile  are  shown  in  Figure  1.  (Tascione,  1988) 

The  Earth’s  chemical  con^sition  spheres,  also  shown  in  Figure  1,  begin  with  the 
homosphere.  This  region  extends  from  the  Earth’s  surface  to  an  altitude  of  approximately 
100  km.  In  the  homosphere,  eddy  diffusion  or  convective  atmospheric  mixing,  keeps  the 
relative  concentrations  of  component  gasses  effectively  constant.  Within  the  homosphere 
lies  the  ozone  layer  at  an  altitude  of  about  30  km.  The  ozone  layer  absorbs  nearly  100%  of 
solar  UV  radiation  between  approximately  2000  A  and  2900  A.  It  is  subsequently  the 
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main  cause  of  atmospheric  warming  in  the  stratosphere.  Above  the  homosphere,  begin¬ 
ning  at  approximately  80-100  km,  is  the  heterosphere.  Here,  the  rate  of  convective  mixing 
drops  off  exponentially  with  altitude.  The  primary  transport  mechanism  in  this  region  is 
therefore  molecular  diffusion.  As  a  result,  the  constituent  species  diffuse  through  the 
region  to  maintain  hydrostatic  equilibrium.  Finally,  the  exosphere  extends  from  approxi¬ 
mately  500  km  to  1000  km  and  above.  In  the  exosphere,  high  energy  neutral  particles  tend 
to  escape  the  Earth’s  atmosphere  due  to  their  high  kinetic  energies  and  the  infrequency  of 
collisions  in  this  region.  The  low  collision  frequency  is  a  consequence  of  the  low  density 
of  atmospheric  particles  and  therefore  the  large  mean  free  path.  (Tascione,  1988) 


Figure  1  Temperature  and  Composition  Spheres  of  the  Atmosphere  (from  Tascione,  1988) 
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B.  MOTIVATION  FOR  STUDY 


1.  Naval  Postgraduate  School  Research 

As  discussed  above,  NPS  has  conducted  two  sounding  rocket  experiments. 
These  experiments  were  made  from  the  White  Sands  Missile  Range  in  New  Mexico  using 
the  MUSTANG  instrument.  They  were  made  in  conjunction  with  the  Naval  Research 
Lab’s  High  Resolution  Airglow  and  Aurora  Spectogr^h  (HIRAAS).  MUSTANG  uses  an 
Ebert -Fastie  spectrograph  to  examine  a  band  from  1800  to  3400  A,  primarily  measuring 
emissions  from  N2  and  NO  with  a  resolution  of  qjproximately  10  A.  The  follow-on  ver¬ 
sion  of  the  MUSTANG  instrument  is  the  higher  resolution  Ionospheric  Spectroscopy  and 
Atmospheric  Composition  (ISAAC)  instrument,  scheduled  for  launch  in  1995  aboard  a 
low  earth  orbit  satellite. 

2.  Motivation  for  High  Resolution  Spectroscopy 

The  10  A  resolution  of  the  MUSTANG  instrument  will  be  greatly  improved  by 
ISAAC  with  a  resolution  of  1.5  A;  however,  there  is  still  a  need  for  new  instruments  with 
even  greater  resolution.  The  All-reflection  Michelson  Interferometer  (AMI)  is  such  an 
instrument.  With  an  expected  resolution  on  the  order  of  a  few  mA,  this  instrument  will 
enable  atmospheric  researchers  not  only  to  obtain  more  accurate  kinetic  temperature  data 
but  also  to  more  closely  study  emission  line  profiles.  In  depth  study  of  emission  line  pro¬ 
files  will  facilitate  a  better  understanding  of  the  physics  of  the  underlying  radiative  trans¬ 
port  mechanisms  (Qeary,  et  al.,  1992).  The  specific  motivation  behind  the  version  of  the 
All-reflection  Michelson  Interferometer  presented  herein  is  the  requirement  to  measure 
the  altitude  distribution  of  atomic  oxygen.  The  AMI  instmment  will  permit  the  direct  mea¬ 
surement  of  the  triplet  profile  at  1304  A.  From  this  measurement,  the  density  of  atomic 
oxygen  in  the  thermosphere  can  be  calculated  without  the  requirement  for  complex  chem¬ 
ical  models  or  absolute  instrument  sensitivity  (Cleary,  et  al.,  1993).  Moreover,  the  design 
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discussed  here  can  be  ad^ted  to  examine  line  profiles  of  numerous  atmospheric  species, 
providing  a  better  understanding  of  the  make-up  of  the  ionosphere.  This  area  of  study  is  of 
particular  interest  to  the  DoD  owing  to  the  importance  of  the  ionosphere  in  the  range  pre¬ 
diction  and  propagation  characteristics  of  over-the-horizon  RADAR  and  long-range  radio 
communications. 
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III.  THEORY 


This  section  introduces  the  basic  theory  required  to  understand  the  method,  results, 
discussion  and  conclusions  regarding  the  All-reflection  Michelson  Interferometer.  A  short 
description  of  the  mathematical  function  used  to  model  the  line  profile  of  interest  is  first 
presented.  A  general  overview  of  interfertnce  and  diffraction  theory  is  then  given  fol¬ 
lowed  by  a  description  of  the  All-reflection  Michelson  Interferometer. 


A.  SOLAR  OXYGEN  1304  A  TRIPLET  LINE  PROFILE 

In  order  to  study  the  characteristics  of  the  AMI  with  respect  to  the  1304  A  triplet,  a 
simple  but  accurate  model  of  that  profile  was  required.  A  straightforward  model  developed 
in  detail  by  Gladstone  (1992)  was  chosen.  This  model  uses  the  sum  of  two  equal  and  offset 
gaussians  to  model  the  self-reversed  spectral  line  profile  according  to  the  equation 


where  nF(x)  is  the  flux  in  ergs  cm'^  sec’^  (km/sec )'^,  x  is  the  Doppler  velocity  in  km/sec 

(converted  to  wavelength  below),  and  tiF  is  the  line  integrated  flux  in 
2  I 

ergs  cm  sec  .  The  quantities  X(,£f  and  x^jj  are  Doppler  velocities  and  are  described  by 
Gladstone  as  a  measure  of  the  offset  of  the  emission  peaks  from  the  center  of  the  self- 
reversed  line  and  a  measure  of  the  dispersion  of  the  line,  respectively. 
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For  the  purposes  of  this  thesis,  x^ff  and  were  converted  from  Doppler  velocity  (in 

km/sec)  to  Doppler  wavelength  shift,  AX  (in  meters),  in  the  following  fashion.  First,  the 
Doppler  wavelength  shift  is  defined  by. 


AX  = 


(2) 


where  x  is  a  Doppler  velocity  as  in  Equation  ( 1 )  and  c  is  the  velocity  of  light.  Xo  is  the  cen¬ 
ter  wavelength  and  X  =  Xq  +  AX.  Therefore, 


X 


_  cAX 


(3) 


Substituting  AX  =  X-X^,  gives 


X  = 


c(X-X^) 


cX 


(4) 


therefore,  the  equations  for  the  offset  and  dispersion  with  respect  to  wavelength  become. 


-c 


(5) 


*dis - Y - 

^0 


(6) 
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Substituting  the  above  two  equations  into  Equation  (1)  gives 


JlF(X)  = 


tcF 


2fn 


/cX 


dis 


V 


-c 


+  e 


(A.-2Xo-h>.o„)' 


.  (7) 


Equation  (7),  giving  flux  as  a  function  of  wavelength,  is  detennined  by  the  parame¬ 
ters  and  given  by  rearranging  Equations  (5)  and  (6)  to  get 


(8) 

(9) 

Equations  (7),  (8)  and  (9)  were  used  to  model  the  emission  line  profile  in  the  IDL  pro¬ 
gramming  language  as  discussed  in  Chapter  V. 


=  Vv+i) 


B.  INTERFERENCE  AND  DIFFRACTION 

Basic  to  an  understanding  of  the  background,  procedures  and  results  of  this  thesis  is  a 
general  understanding  of  the  concepts  of  interference  and  diffraction  and  their  tq^plicabU- 
ity  in  the  area  of  spectroscopy.  Interference  and  diffraction  theory  are  presented  separately 
below.  More  complete  coverage  of  these  two  subjects  is  available  in  any  number  of  stan¬ 
dard  texts  on  basic  physics,  electromagnetic  theory,  or  optics  such  as  Pedrotti  and  Pedrotti 
(1987). 
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1.  Interference  Theory 

Interference  can  be  described  as  the  enhancement  or  diminution  of  the  intensity 
of  electromagnetic  radiation  caused  by  the  superposition  of  two  or  more  individual  waves. 
The  classical  interference  of  two  waves  causes  alternate  enhancement  and  reduction  of 
wave  intensity  in  what  can  be  described  as  an  interference  or  fringe  pattern. 

In  classical  two  beam  interference,  two  electromagnetic  waves,  represented  by 
their  vector  components  Ej  and  emanate  from  a  single  source,  travel  along  different 
paths  and  then  recombine  at  some  point  P  identified  by  its  position  vector  r.  Assuming  that 
the  waves  are  of  equal  amplitude,  they  can  be  described  by 


E|  =  Eq  cos  (k|  •  r  -  ojt) 

(10) 

E2  =  Eo  cos  (k2  •  r  -  cot  + 

(11) 

where  Eq  represents  the  amplitude;  the  initial  phase  difference;  (O,  the  wave  frequency; 
and  k,  the  propagation  vector.  At  the  point  of  recombination,  P,  the  electromagnetic  field  is 
given  by  the  superposition  of  Ej  and  E2, 

Ep=E|  +  E2.  (12) 

In  optics,  electromagnetic  waves  are  mote  conveniently  represented  by  their  irradiance  I, 
the  value  that  is  commonly  measured.  Irradiance  is  related  to  the  electric  field  vector  by 
I  =  Cqc  (E^),  with  Cq  the  permittivity  constant,  c  the  speed  of  propagation,  and  the  angled 
brackets  indicating  the  time  average  of  the  quantity  inside.  At  point  P  the  irradiance  is. 

Ip  =  Cqc  (Ep^)  =  Efl  c  (Ej^  +  +  2Ei  •  E2),  (13) 
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or, 


Ip  "  I|  +  I2  +  Ii2.  (I^) 

where  Ix  and  I2  are  the  irradiances  of  the  constituent  waves  and  1x2  ^  ^  additional  irradi- 
ance  term  due  to  the  interference.  This  interference  term  is  then 

Ix2  =  2eoc<Ex‘E2>.  (15) 

Through  the  use  of  trigonometric  identities  and  the  time  average  characteristics  of  the  trig¬ 
onometric  functions  below. 


(  sin^  (cat) )  =  1/2  (16) 

<cos2((ijt)>=l/2  (17) 

(  sin  (<nt)  cos  (cat) )  =  0,  (18) 

the  time  average  term  of  the  interference  becomes 

(  Ex  •  E2  >  =  V2  Eo^  cos  [(kx  -  k2)  •  r  +  <|>],  (19) 


where  the  term  in  square  brackets  represents  the  phase  difference  between  the  constituent 
waves,  hereafter  t  noted  by  5.  The  interference  term  of  the  irradiance  is  then 

1x2  =  Eo  c  Eq^  cos  5.  (20) 

Using  the  time  average  characteristic  of  cos^(cat)  above,  the  irradiance  terms  representing 
the  original  two  waves  become 


(21) 
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Therefore  the  interference  term  becomes. 


Il2  =  2 1  cos  5  (22) 

and  the  irradiance  at  point  P  is  then 

Ip  =  2 1  +  2 1  cos  6.  (23) 

The  irradiance  therefore  becomes  a  maximum  when  cos  5  is  +1  or  5  =  2nm,  (m  =  0,  ±1 , 

±  2, ...)  and  a  minimum  when  cos  d  is  -1  or  5  =  (2m-f-l)  n,  (m  =  0,  ±1,  ±  2, ...).  (Pedrotti 
and  Pedrotti,  1987) 

As  described  above,  two  beam  interference  is  the  periodic  enhancement  and 
reduction  of  the  irradiance  of  two  superimposed  waves.  It  is  caused  by  the  phase  shift,  5, 
between  the  waves  and  can  occur  when  the  planes  of  constant  phase  (or  wavefronts)  are 
split  among  different  paths.  The  total  phase  shift,  6,  is  due  to  the  combination  of  the  initial 
phase  shift  of  the  waves  and  the  difference  in  path  length  between  the  two  beams.  Two 
beam  interference  theory  can  be  generalized  to  describe  wavefront  division  interference  of 
multiple  waves. 

A  final  important  aspect  of  interference  is  the  concept  of  coherence  length.  Opti¬ 
cal  waves  do  not  travel  in  perfect  and  infinite  sinusoids.  They  travel  in  groups  that  appear 
sinusoidal  for  a  finite  length  and  then  discontinuously  change  phase^as  the  next  group 
begins.  The  average  lifetime  of  the  group  is  called  the  coherence  time;  the  average  length 
of  this  coherent  pulse  is  refened  to  as  the  coherence  length  and  varies  with  the  source  of 
electromagnetic  energy  (Pedrotti  and  Pedrotti,  1987).  Its  importance  lies  in  the  fact  that 
for  interference  to  take  place,  paired  groups  of  waves  must  arrive  together  at  the  field 
point,  P.  If  the  path  length  difference  is  greater  than  the  coherence  length,  paired  groups  do 
not  arrive  within  the  coherence  time,  and  interference  does  not  occur.  Throughout  the  con- 
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duct  of  this  thesis,  coherence  length  was  arbitrarily  and  artificially  set  to  a  very  high  value 
such  that  it  was  effectively  infinite. 

2.  Diffraction  Theory 

Diffraction  is  similarly  an  effect  of  the  wave  nature  of  light.  Diffraction  is 
described  by  Pedrotti  and  Pedrotti  (1987)  as  “...any  deviation  from  geometrical  optics  that 
results  from  the  obstruction  of  a  wavefront  of  light.”  Interference  between  diffracted 
waves  produces  a  pattern  of  fringes.  As  shown  in  Figure  2  below,  the  diffraction  of  light 
by  multiple  slits  produces  a  fringe  pattern  consisting  of  bright,  distinct  and  widely  sepa¬ 
rated  maxima.  A  diffraction  grating  is  a  device  containing  evenly  spaced,  multiple  slits  or 
grooves  that  takes  advantage  of  the  sensitivity  of  its  diffraction  pattern  to  the  wavelength 
of  the  incident  light. 


Figure  2  Pattern  for  Multiple  Slit  Fraunhofer  Diffraction  (Pedrotti  and  Pedrotti,  1987) 


Far  field  or  Fraunhofer  diffraction  is  diffraction  where  the  observation  plane  is  a 
large  distance  from  the  diffraction  source  such  that  the  wavefronts  can  be  considered  pla- 
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nar.  In  classic  Fraunhofer  diffraction,  multiple  slit  diffraction  using  a  grating  produces 
irradiance  maxima  (as  shown  in  Figure  2  above)  according  to  the  grating  equation 

mX  =  a  ( sin0.  +  sin9^) .  (24) 

Here  m  is  the  order  of  diffraction,  0|  and  Og,  are  the  incident  and  diffraction  angles  mea¬ 
sured  from  the  grating  normal,  and  a  is  the  distance  between  slits  on  the  grating.  Usually 
the  grating  ruling  density,  p  =  1/a,  is  given  and  the  grating  equation  then  becon^s 

pmX  =  ( sinO.  +  sin9^) .  (25) 

Therefore,  a  diffraction  grating  takes  a  beam  of  monochromatic  incident  light  and,  upon 
reflection  from  the  grating,  splits  it  along  multqile  paths  according  to  the  grating  equation. 
The  zero  order  diffraction  path  is  the  same  for  all  wavelengths  and  is  merely  a  reflection 
following  Snell’s  law.  Positive  order  diffractions  diffract  light  at  an  angle  greater  than  that 
of  the  zero  order  and  negative  order  diffractions  di^act  light  at  an  angle  less  than  that  of 
the  zero  order.  As  can  be  seen  in  Figure  2  above,  the  intensity  of  the  diffracted  beam  falls 
off  as  the  diffraction  angle  (and  thus  diffraction  order)  increases  in  absolute  value. 

Figure  3  shows  an  illustration  of  a  diffraction  grating  with  multiple  orders  of  diffraction. 
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Figure  3  Reflection  Grating  with  Multiple  Diffraction  Orders 


There  is  no  fundamental  physical  difference  between  interference  and  diffrac¬ 
tion  patterns.  Fringes  caused  by  the  diffraction  of  light  are  more  rigorously  termed  inter¬ 
ference  fringes,  but  are  commonly  called  diffraction  patterns.  The  words  interference  and 
diffraction  are  often  used  interchangeably;  however,  in  this  thesis,  interference  will  refer 
to  the  effect  producing  a  fringe  pattern  while  diffraction  will  refer  to  the  change  in  direc¬ 
tion  of  a  beam  caused  by  an  obstruction  or  aperture. 


C.  THE  ALL-REFLECTION  MICHELSON  INTERFEROMETER 

An  interferometer  is  a  device  that  splits  and  then  recombines  light  to  create  a  pattern 
of  interference  fringes  for  a  specific  purpose.  Generally  classified  by  the  method  used  to 
separate  the  incident  beam  of  light,  interferometers  are  termed  either  an^litude  or  wave- 
front  splining  devices.  This  section  discusses  the  basic  theory  behind  the  plane  grating 
All-reflection  Michelson  Interferometer.  The  AMI  is  classed  as  an  amplitude  splitting 


15 


device  because  of  the  use  of  a  diffraction  grating  as  a  beam  splitter  to  divide  the  input 
beam  among  separate  paths. 

The  general  design  for  the  All-reflection  Michelson  Interferometer  examined  in  this 
thesis  was  developed  in  detail  by  Wallace  (1992).  This  basic  all-reflection  interferometer 
design  consists  of  an  aperture,  a  plane  diffraction  grating,  two  plane  mirrors,  and  a  detec¬ 
tor.  The  grating  splits  the  input  light  into  two  paths  corresponding  to  the  zero-  and  minus- 
one-order  diffraction.  Each  path  travels  to  a  plane  mirror  which  reflects  the  light  back  to 
the  grating.  The  light  again  reflects  off  the  grating  utilizing  the  minus-one  and  zero  diffrac¬ 
tion  orders,  respectively.  The  light  finally  recombines  along  a  commm  path  from  the  grat¬ 
ing  to  the  detector.  The  input  light  is  assumed  to  be  monochromatic  and  collimated  for  the 
purposes  of  this  discussion.  The  general  layout  of  the  elements  of  the  All-reflection  Mich¬ 
elson  Interferometer  is  shown  in  Figure  4.  Figure  5  shows  the  path  of  several  rays  in  die 
AMI. 
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Figure  5  Ray  Paths  in  a  Typical  All-reflection  Michelson  Interferometer. 


In  the  AMI,  the  path  length  measured  along  each  of  the  two  paths  (hencefoith  known 
as  the  “zero-order  path”  and  the  “minus-one-order  path”)  is  equal.  When  the  wavelength 
of  the  input  light  is  equal  to  the  “tuned”  wavelength  (known  as  Xq)*  *he  wavefronts  of  the 
two  recombined  beams  are  parallel.  The  interference  at  the  detector  is  constructive  at  all 
points  and  the  interference  pattern  has  zero  spatial  frequency.  When  the  wavelength  of  the 
incident  light  is  offset  from  Ao,  however,  a  shift  in  the  diffraction  angles  from  that  of  the 
tuned  wavelength  results.  This  shift  causes  a  tilt  between  the  wavefronts  of  the  two  beams. 
The  spatial  frequency  of  the  resulting  interference  pattern  is  wavelength  dependent.  The 
fourier  transform  of  this  spatial  fringe  pattern  indicates  the  amount  of  wavelength  offset, 
AX,  from  the  tuned  wavelength,  where  AX  =  X  -  Xq.  The  instrument  presented  here,  how¬ 
ever,  is  incapable  of  differentiating  between  equal  offsets  both  greater  and  less  than  Xq.  A 
more  detailed  treatment  of  the  so-called  spatial  heterodyning  effect  can  be  found  in  Wal¬ 
lace  (1992). 
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Recall  that  the  AMI  is  designed  so  that  the  length  of  each  of  the  two  paths  is  equal. 
The  angles  of  incidence  and  diffraction,  and  the  placement  of  the  mirrors  are  detennined 
as  described  below.  In  the  AMI  design,  use  of  a  plane  difu'action  grating  serves  to  diffract 
the  incident  light  along  paths  according  to  the  general  grating  equation  (Equation  (25)), 
repeated  here. 


pmA,  =  sin0. +  sin0^,  (26) 

where  again  p  is  the  grating  density  (lines  Anm),  m  is  the  diffraction  order,  A.  is  the  wave¬ 
length  (mm),  and  0j  and  0q,  are  die  incident  angle  and  the  m*^  order  diffraction  angle, 
respectively.  All  angles  are  measured  clockwise  from  the  grating  normal.  In  order  that  the 
zero-  and  minus-one-order  paths  have  equal  lengths,  their  respective  mirrors  must  be  equi¬ 
distant  from  the  grating. 

Multiple  diffraction  orders  emanating  from  the  grating  reduce  the  overall  intensity  of 
light  incident  upon  the  detector.  As  the  instrument  only  uses  the  zero-  and  minus-one- 
order  initial  diffractions,  the  light  incident  upon  the  detector  can  be  maximized  by  elimi¬ 
nating  all  diffractions  of  higher  order.  This  is  done  by  proper  selection  of  p  with  respect  to 
A,  so  that  0}  and  all  higher  order  diffraction  angles  are  greater  than  or  equal  to  90°  and  0.2 
and  all  lower  order  diffraction  angles  are  less  than  or  equal  to  -90°  (Cleary,  et  al.,1992). 

Setting  0|  equal  to  90°  and  substituting  into  Equation  (26)  above  gives, 

pA,  =  sinO.  +  1  (27) 

or, 

0j  =  asin  (pA,-  1) .  (28) 
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Similarly,  setting  O.2  to  -90^  and  substituting  into  the  grating  equation  gives 

-2pA,  =  sin0.  -  1  (29) 

or, 

Gj  =  asin(-2pX+l) .  (30) 

Substituting  Equation  (28)  iiito  Equation  (30)  and  solving  for  p  gives 

p  =  ^.  (31) 

Equation  (31)  gives  the  optimum  ruling  density  with  respect  to  wavelength  for  mini¬ 
mum  loss  of  light  intensity  due  to  wasted  diffraction  orders.  In  practice,  one  selects  the 
closest  commercially  available  grating  with  a  grating  density  greater  than  or  equal  to  that 
given  in  Equation  (31)  above.  (Wallace,  1992) 

In  fact,  as  pointed  out  by  Wallace  (1992),  selection  of  a  grating  density  precisely 
equal  to  that  given  above  would  force  the  incident  and  minus-hrst-order  diffraction  angles 
to  be  identical,  an  undesirable  characteristic  for  this  q>plication  as  the  minus-one-order 
mirror  would  then  be  in  the  path  between  the  grating  and  source.  Therefore,  Equation  (31) 
has  been  established  by  Wallace  (1993)  as  the  minimum  acceptable  grating  density  for  the 
elimination  of  wasted  orders  at  the  first  pair  of  diffractions.  Once  selection  of  the  grating 
density,  wavelength  of  interest,  and  incident  angle  has  been  completed,  the  positions  of  the 
components  of  the  interferometer  are  determined  by  the  grating  equation.  For  the  zero 
order  mirror. 


(32) 
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and  for  the  minus  first  order  mirror. 


0_j  =  asin  (-pA.- sin0.) .  (33 

The  final  position  of  each  of  the  components  can  be  made  after  choosing  tf^  distances 
from  the  grating  to  the  source,  mirrors,  and  detector.  This  selection  can  be  based  almost 
entirely  on  size  constraints  of  the  completed  instrument.  Selection  of  die  angle  of  inci¬ 
dence  will  be  treated  in  Chqiter  V  of  this  thesis. 


IV.  RAY  TRACING  PROCEDURES 


The  method  of  ray  tracing  was  used  to  model  interference  patterns  for  several  ver¬ 
sions  of  the  proposed  instrument.  Additionally,  once  the  instrument  layout  was  finalized, 
numerous  interference  simulations  were  performed  to  determine  some  final  performance 
parameters  of  the  instrument.  This  chs^ter  contains  a  short  description  of  the  software 
used  for  ray  tracing  and  interference  pattern  constraction  and  describes  limitations  of  the 
software  for  the  present  application. 

A.  DART  RAY  TRACING  PROGRAM 

The  majority  of  the  results  presented  in  this  thesis  were  developed  by  simulating  the 
instrument  using  the  DART  ray  tracing  program  of  Atkinson  (1993).  DART  is  a  Windows- 
based  PASCAL  program  written  for  the  IBM  PC.  It  was  specifically  designed  for  complex 
optical  ray  tracing  t^lications  and  is  described  in  detail  by  Atkinson  (1993).  In  DART, 
instiuments  are  modeled  by  combining  components  from  an  extensive  element  library. 
Apertures  in  the  DART  program  (hereafter  called  DART  {qrertures)  are  modeled  using  a 
multi-function  user  interface  and  can  be  tailored  in  both  shape  and  size  to  the  specific 
application.  The  light  source  can  be  monochromatic  radiation  of  any  wavelength  or  can  be 
composed  of  multiple  wavelengths  at  selectable  intensities. 

DART  uses  two  three-dimensional  right-hand  coordinate  systems  for  the  placement 
of  optical  elements.  The  first,  called  the  “global”  coordinate  system  is  used  for  the  place¬ 
ment  of  instrument  elements  in  relation  to  the  instrument  aperture  which  is  usually  at  the 
origin.  The  initial  direction  of  the  input  beam  is  in  the  positive  Z  direction.  The  most  com¬ 
mon  layout  of  a  planar  instrument  is  to  place  ail  elements  in  the  X-Z  plane.  The  second 
coordinate  system  is  the  “local”  coordinate  system.  This  set  of  coordinates  is  fixed  with 
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respect  to  each  element  in  the  system.  The  origin  of  each  local  coordinate  system  is  cen¬ 
tered  on  the  element  with  the  z-axis  in  the  direction  of  the  outward  normal  of  the  element 
surface.  Translation  and  rotation  of  an  element  changes  its  position  and  direction  in  the 
global  coordinate  system  but  not  in  the  local  system.  Throughout  this  thesis,  axes  and 
directions  described  by  ci^ital  letters  will  denote  global  coordinates  while  lower  case  let¬ 
ters  will  indicate  the  local  coordinates.  An  illustrative  two-element  instrument  with  the 
coordinate  systems  indicated  is  shown  in  Figure  6. 

The  DART  program  produces  a  rough  layout  of  the  instnunent  being  modeled  which 
can  be  viewed  from  any  angle.  If  desired,  DART  traces  the  ray  paths  on  this  layout  dis¬ 
play.  DART  layout  drawings  only  trace  those  rays  that  contact  the  next  element  in 
sequence — missed  rays  are  not  drawn.  Plane  surfaces  may  be  used  to  mask  off  undesired 
rays  anywhere  in  the  instrument.  The  final  plane  surface  in  any  instrument  serves  as  the 
instrument’s  detector.  The  specific  layout  for  the  AMI  is  described  in  detail  in  Ch£q>ter  V. 

The  interference  patterns,  or  interferograms,  produced  by  DART  are  viewed  in  an 
interference  “window”  on  the  computer  screen.  An  interference  profile,  created  by  verti¬ 
cally  summing  the  400  element  wide  interferogram,  may  be  saved  as  a  binary  data  file. 
The  result  is  a  file  consisting  of  400  pairs  of  intensity-versus-detector-x -position  values. 
Additionally,  an  instrument  parameter  file  may  be  created  which  contains  a  complete  list¬ 
ing  of  the  input  parameters  for  the  selected  instrument  path. 
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Figure  6  Global  and  Local  Coordinate  Systems 


B.  RAY  TRACING  SOFTWARE  LIMITATIONS 

Although  of  inestimable  benefit,  the  DART  program  contains  several  features  which 
limit  its  usefulness  in  the  present  application.  The  first  of  these  is  the  limit  of  400  bins  in 
which  to  sum  the  interference  data.  This  400  bin  limit  constrains  the  resolution  demonstra¬ 
ble  by  DART  interference  files  compared  to  a  512  or  1024  pixel  one.-dimensional  array 
detector  (two  commonly  available  sizes).  The  number  of  bins  in  DART  was  based  on  the 
limit  of  screen  resolution  on  a  typical  computer  monitor.  A  suggested  in^rovement  is  to 
include  an  option  to  select  greater  numbers  of  bins  and  have  DART  calculate  ray  paths  and 
interference  patterns  at  a  higher  resolution  than  can  be  displayed  on  the  screen.  As  most  of 
the  analysis  of  ray  tracing  results  will,  in  all  likelihood,  be  done  using  other,  more  robust 
calculating  programs  such  as  IDL  or  FORTRAN,  this  change  would  be  beneficial  to  most 
DART  users. 
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A  second  limitation  of  the  DART  program  is  the  requirement  that  the  number  of  rays 
striking  the  detector  be  identical  in  the  two  instrument  paths  (for  multiple  wavelength 
interference).  If  the  number  of  rays  at  each  path’s  detector  is  unequal,  interference  is  not 
calculated.  Although  not  troublesome  when  using  a  small  aperture  and  beam  size  com¬ 
bined  with  large  components  and  a  large  detector,  the  limitation  is  often  iqjparent  when 
using  small  components.  A  more  beneficial  approach  would  have  been  for  DART  to  dis¬ 
count  rays  from  one  path  that  can’t  be  matched  with  rays  from  the  other,  while  allowing 
the  interference  to  be  computed.  DART  could  then  aleit  the  user  to  the  number  of  unusable 
rays. 

The  final  limitation  of  note  is  the  requirement  that  rays  emanating  from  the  instru¬ 
ment  ^rture  initially  travel  in  the  positive  Z  direction.  This  requirement  can  sometimes 
force  the  instrument  designer  to  place  a  majority  of  critical  components  on  the  negative  Z 
side  of  the  X-axis.  Although  easily  woiked  around,  this  complicates  the  translation  from 
design  on  DART  to  placement  of  actual  elements  on  an  optical  bench. 
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V.  EXPERIMENTAL  PROCEDURE  AND  RESULTS 


A.  IDL  FOURIER  TRANSFORM  PROCEDURES 

In  the  AMI,  the  spatial  frequency  of  the  interference  pattern  is  wavelength  dependent. 
One  recovers  the  spectrum  of  the  incident  light  source  by  taking  the  fourier  transform  of 
the  interference  profile.  Interactive  Data  Language  (IDL)  for  PC  with  Windows  was  used 
in  order  to  examine  the  spectrum  of  die  fringe  patterns  created  by  the  DART  af^lication. 
An  IDL  procedure  was  written  to  read  interference  data  from  a  DART  output  file  and  dien 
take  the  fast  fourier  transform  of  that  data  using  the  IDL  FFT  function.  The  procedure 
plots  the  interference  profile  (relative  intensity  versus  the  x-axis  position  of  the  detector) 
and  the  fast  fourier  transform  spectrum  (relative  intensity  versus  position  in  the  FFT 
array).  IDL  places  the  fourier  transformed  data  into  an  array  of  the  same  dimensions  as  the 
input  data  and  then  plots  the  full  FFT  spectrum  consisting  of  both  positive  and  negative 
(ghost)  frequencies. 

For  the  purposes  of  this  thesis,  the  full  FFT  spectrum  contains  redundant  information. 
A  so-called  wavelength  spectrum  was  obtained  by  taking  one  half  of  the  full  FFT  spec¬ 
trum  and  converting  the  horizontal  axis  to  wavelength.  For  clarity  of  presentation,  positive 
frequencies  of  the  fourier  transform  were  used  for  incident  wavelengths  greater  than  Xq 
and  negative  frequencies  for  wavelengths  less  than  X^.  In  the  present  case,  the  fourier 
spectrum  consists  of  400  elements,  half  of  -viiich  are  used  to  plot  the  wavelength  spec¬ 
trum.  FFT  spectra  were  converted  to  wavelength  spectra  by  another  IDL  procedure.  These 
IDL  plotting  procedures  and  functions  are  also  presented  in  Appendix  A. 
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B.  EXPERIMENT  DOCUMENTATION 


The  research  conducted  as  part  of  this  thesis  produced  three  distinct  foiins  of  output. 
The  DART  interferogram  is  a  two-dimensional  representation  of  the  interference  pattern 
at  the  instrument  detector.  Each  ray  striking  the  detector  is  plotted  using  a  gray  scale  based 
on  the  intensity  of  the  interfering  light.  The  interferogram  is  a  representation  of  how  inter¬ 
ference  would  appeal  were  the  light  source  visible.  The  second  form  of  output  is  the  inter¬ 
ference  profile.  This  is  a  plot  of  intensity  versus  x-axis  position;  it  represents  the  pattern 
that  would  be  obtained  from  a  one-dimensional  array  detector.  The  final  and  most  useful 
form  of  output  is  the  wavelengdi  spectrum  (or  alternatively  the  full  FFT  spectrum).  Once 
this  plot  is  calibrated  for  wavelength,  the  peaks  directly  translate  to  the  difference,  AX, 
between  input  wavelength  and  tuned  wavelength. 

The  quality  of  each  of  the  plots  listed  above  can  be  enhanced  simply  by  increasing 
the  number  of  rays  produced  at  the  DART  :q)erture.  However,  speed  of  computation 
within  DART  is  inversely  related  to  the  number  of  rays  produced;  therefore,  most  interfer¬ 
ence  runs  were  completed  using  an  intermediate  number  of  rays.  Documentation  for  this 
thesis  was  later  prepared  by  re-running  desired  simulations  using  a  greater  number  of  rays. 

In  DART,  intensities  are  relative  values  determined  by  the  number  of  rays  striking  the 
detector.  These  intensities  can  be  scaled  by  whole  numbers  when  entering  multiple  wave¬ 
lengths  but  absolute  intensities  are  difficult  to  introduce.  Additionally,  intensities  for  line 
profiles  calculated  using  Equation  (7)  have  been  scaled  and  rounded-off  for  ease  of  entry 
into  DART.  Since  relative  intensities  are  sufficient  for  this  thesis,  no  attempt  was  made  to 
introduce  absolute  intensities  herein.  Comparison  of  relative  intensities,  however,  was 
used  extensively  in  examining  the  behavior  of  the  AMI. 
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C.  RAY  TRACING  PROGRAM  SETUP 

The  setup  of  the  AMI  instniment  in  the  ray  tracing  program  was  undertaken  in  three 
steps.  First,  the  interferometer  portion  alone  was  modeled  using  a  rectangular,  collimated 
beam  of  monochromatic  light  as  the  source  at  the  DART  program  q}erture.  Then  an 
instrument  aperture  more  closely  resembling  the  planned  AMI  design  was  modeled  using 
an  off-axis  partdwlic  reflector  with  the  DART  tq>erture  at  the  focus  of  the  paraboloid. 
Finally,  these  two  DART  qiplications  were  combined  to  form  a  model  of  the  AMI.  Addi¬ 
tional  elements  are  expected  to  be  included  in  later  models  of  the  actual  flight-qualified 
instrument  in  order  to  better  collect  and  direct  the  light  into  the  spectrograph  ^rture. 

Once  the  complete  instrument  was  modeled,  several  different  versions  of  the  DART 
aperture  were  created.  DART’s  aperture  determines  the  form  of  the  light  incident  on  the 
instrument’s  first  element.  The  first  version  contained  an  ideal  point  source  which  con¬ 
sisted  of  rays  emanating  from  a  single  point  and  spreading  spherically  from  the  DART 
aperture  to  the  surface  of  the  parabolic  reflector.  The  second  version  used  a  circular  “pin¬ 
hole”  source  ranging  in  diameter  from  100  pm  to  500  pm.  The  final  version  consisted  of  a 
vertical  slit  aperture  of  varying  widths  and  a  height  of  one  centimeter.  In  both  the  pinhole 
and  the  slit,  rays  emanated  and  spread  spherically  outward  from  the  entire  surface  area  of 
the  source.  A  more  complete  description  of  the  three  instrument  setups  is  provided  below. 

1.  Interferometer 

As  described  above,  the  interferometer  section  of  the  AMI  instrument  was  mod¬ 
eled  with  the  DART  ray  tracing  application.  The  model  consisted  of  a  square  {q)erture,  a 
square  diffraction  grating,  two  plane  circular  mirrors,  and  a  plane  surface  for  use  as  a 
detector.  The  overall  dimensions  of  the  instrument  were  selected  to  be  representative  of 
those  in  a  rocket  or  space-based  version.  The  distances  from  the  diffraction  grating  to  the 
various  elements  are  shown  in  Table  I. 
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TABLE  I.  INTERFEROMETER  DIMENSIONS 


Element 

Distance  to 
Grating 

Apertiue 

40  cm 

Mirrors 

20  cm 

Detector 

25  cm 

The  grating  density  of  the  instrument  was  then  selected  by  using  Equation  (31) 
to  determine  the  optimal  minimum  grating  density.  Substituting  in  the  desired  tuned  wave¬ 
length  of  1304  A  gives 


P  = 


A  -  2 

3  (1.304xl0~^mm) 


=  5112 


lines 
mm  ’ 


(34) 


which  leads  to  a  commonly  available  grating  choice  of  p  =  6000  lines/mm.  The  initial  size 
of  each  of  the  elements  in  the  DART  simulation  was  chosen  to  be  large  enough  so  that  all 
light  from  the  DART  s^rture  would  strike  each  element  in  turn  and  finally  strike  the 
detector.  In  that  way,  accuracy  in  placement  and  orientation  of  each  element  was  less  crit¬ 
ical.  The  elements  were  placed  and  oriented  according  to  the  theory  described  earlier 
using  the  dimensions  specified  in  Table  I  and  the  selected  incident  angle.  Calculations  and 
initial  conversion  into  DART’s  coordinate  system  were  simplified  by  using  a  MATHCAD 
document  written  by  Wallace  (1992).  Once  entered  into  the  DART  program,  interference 
patterns  were  created  by  running  DART’s  interference  macro  for  a  number  of  different 
wavelengths  of  monochromatic  light.  Recall  that  the  DART  program  requires  separate 
instrument  layouts  for  each  path.  An  isometric  view  of  each  path  of  the  interferometer 
with  one  ray  displayed  is  shown  in  Figure  7  and  8.  The  resulting  interferogram  for  an 
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input  wavelength  of  1304.02  A  is  shown  in  Figure  9.  All  three  figures  were  drawn  by  the 
DART  program.  Appendix  B  includes  the  DART  parameter  file  for  each  path  shown. 


Figure  7  Interferometer  Layout  of  the  0  Order  Path  Produced  by  DART  (not  to  scale) 
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2.  OfT-axis  Parabolic  Reflector 


The  off-axis  parabolic  reflector  was  used  as  a  collimating  element  between  the 
source  and  the  grating.  The  paraboloid  was  modeled  in  DART  according  to  the  specifica¬ 
tions  for  the  Oriel  model  45347  parabolic  reflector.  This  particular  reflector  was  used  by 
Risley  (1993)  in  determining  the  physical  problems  associated  with  constructing  the  AMI 
for  use  in  the  visible  and  near  ultraviolet  spectrum. 

In  order  to  create  an  off-axis  parabola  in  the  DART  iqrplication,  the  full  para¬ 
bolic  mirror  must  first  be  entered.  Then  a  plane  surface  mask  must  be  placed  between  the 
parabola  and  the  next  element  in  sequence.  The  purpose  of  the  mask,  which  is  transparent 
to  all  rays  striking  it,  is  to  eliminate  all  but  those  rays  which  would  reflect  off  a  real  off- 
axis  parabola.  This  is  accomplished  because  DART  only  traces  rays  that  strike  the  next 
element  in  sequence  (see  Chapter  IV). 

In  order  to  accurately  construct  the  DART  version  of  the  off-axis  parabola,  its 
orientation  must  be  such  that  no  rays  are  required  in  the  negative  Z  direction  of  DART’s 
global  coordinate  system  (see  Ray  Tracing  Software  Limitations  above).  Therefore,  the 
parabola  was  oriented  as  shown  in  Figure  10  below.  This  figure  shows  the  general  layout 
of  the  off-axis  parabola  and  includes  the  trace  of  a  small  number  of  rays. 

As  can  be  seen  in  Figure  12  and  in  several  of  the  interferograms,  interference 
profiles  and  wavelength  spectra  that  foUow,  the  off-axis  parabola  introduces  an  increase  in 
relative  intensity  of  light  on  one  side  of  the  detector.  This  rise  in  intensity  is  due  to  a  natu¬ 
ral  increase  in  the  number  of  rays  leaving  the  parabola  mask  on  the  side  farthest  away 
from  the  focus.  Were  the  absolute  intensity  significant  in  this  thesis,  the  rise  would  have  to 
have  been  accounted  for  in  subsequent  analysis. 
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Figure  10  Off-axis  Parabolic  Reflector  Simulation 


3.  All-Reflection  Michelsun  Interferometer 

The  two  sections  of  the  instrument  described  above  were  combined  into  a  single 
model  forming  the  current  version  of  the  AU-reflection  Michelson  Interferometer.  This 
instrument  was  initially  designed  with  all  elements  sized  larger  than  required  to  ensure 
that  all  rays  would  strike  the  detector.  After  the  model  was  rotated  for  the  addition  of  the 
off-axis  parabola,  the  interferometer  section  was  placed  20  cm  further  from  the  DART 
aperture  to  facilitate  adding  additional  optics  to  the  simulation  if  required.  The  size  of  each 
element  in  the  instument  was  then  reduced  and  the  shape  changed  to  more  closely  match 
commonly  available  components  or  components  planned  for  use  in  other  current  NPS  ion¬ 
ospheric  research  projects  (Risley,  1993).  The  size  and  shape  of  the  chosen  components 
are  shown  in  Table  II. 
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TABLE  n.  COMPONENTS  USED  IN  DART  AMI  SIMULATION 


Component 

Shape 

Size  (cm) 

Grating 

Square 

5.0  X  5.0 

Mirrors 

Circular 

d  =  2.50 

Detector 

Rectangular 

2.50  X  1 .00 

The  complete  instrument  consists  of  two  DART  parameter  files,  each  represent¬ 
ing  one  path  taken  by  the  rays.  These  hies  are  also  included  in  Appendix  B.  The  complete 
layout  of  the  AMI  used  for  all  further  DART  interference  runs  is  shown  in  Figures  1 1 
and  12.  Sample  ray  traces  are  also  shown. 


Figure  1 1  AMI  DART  Layout  and  Ray  Trace  of  the  Zero-Order  Path 
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Figure  12  AMI  DART  Layout  and  Ray  Trace  of  the  Minus-One-Order  Path 


D.  AMI  MODEL  VALIDATION  AND  OPTIMIZATION 

The  steps  taken  to  validate  the  AMI  design  were  divided  into  five  major  sets  of 
DART  interference  runs.  The  first  set  was  performed  using  an  ideal  point  source  aperture. 
Runs  were  completed  using  both  monochromatic  and  multiple  wavelength  sources.  The 
aperture  was  then  altered  by  offsetting  the  source  from  the  focus  of  the  ofF-aais  parabola. 
A  “real  point  source”  was  then  used  in  the  simulation  in  order  to  note  the  .jffects  of  using  a 
realistic  optical  pinhole  at  the  focus  of  the  off-axis  parabola.  This  is  the  element  most 
likely  to  be  used  in  the  actual  construction  of  the  AMI  (Risley,  1993).  The  pinhole  was 
then  elongated  to  form  a  slit,  another  viable  option  for  the  actual  instrument.  Finally,  the 
ideal  point  source  was  again  used  for  runs  conducted  to  simulate  the  input  of  an  emission 
line  profile. 
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1.  Ideal  Point  Source  with  Monochromatic  Incident  Light 

An  ideal  point  source  was  entered  as  the  aperture  of  the  DART  parameter  file  by 
setting  all  aperture  dimensions  equal  to  zero  and  selecting  an  opening  one  point  across  in 
both  the  X  and  Y  directions.  It  was  first  verified  that  an  input  wavelength  equal  to  the 
tuned  wavelength  (X 1=  =  1304  A)  produced  no  fringe  pattern.  Although  there  were  no 

fringes  visible  on  the  interferogram,  the  wavelength  domain  was  examined  to  ensure  that 
no  peaks  were  present.  That  was  indeed  the  case.  Numerous  values  of  AX.  were  then 
entered  and  the  results  examined.  Of  specific  ifiterest  were  the  minimum  and  maximiun 
allowed  values  of  AX.  which  could  be  determined  through  analysis  of  the  fourier  trans¬ 
formed  interference  patterns. 

The  first  example  demonstrates  a  typical  interference  pattern  and  profile  along 
with  the  resulting  fourier  spectrum.  Figures  x3  through  IS  show  the  resulting  interfero¬ 
gram,  interference  profile,  and  FFT  spectrum,  respectively,  for  an  incident  wavelength  of 
1304.015  A.  The  minima  in  the  interference  profile  are  at  zero  and  therefore  the  profile  is 
not  centered  about  the  x-axis.  This  produces  a  large  DC  term  at  X^  in  the  fourier  plot, 
which  is  present  in  Figure  15.  This  interference  run  was  conducted  using  a  block  of 400  by 
400  rays  sent  from  the  aperture  toward  the  off-axis  parabola.  Figures  16  through  18  show 
the  same  data  using  a  120  by  120  ray  aperture  to  demonstrate  the  decrease  in  the  “signal- 
to-noise”  ratio  that  results  from  using  a  smaller  number  of  rays.  All  ^bsequent  ideal  point 
source  plots  use  a  400  x  400  ray  aperture  unless  otherwise  labeled. 


t 
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Relative  intensity 


Figure  17  Ideal  Point  Source  Interference  Profile  for  1304.01 5  A  (120  x  120  rays) 


Figure  18  Ideal  Point  Source  FFT  Spectrum  for  1304.015  A  (120  x  120  rays) 
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The  next  series  of  results,  shown  in  Figures  19  through  2 1 ,  demonstrates  inter¬ 
ference  using  a  value  of  AA.  near  the  center  of  the  expected  bandwidth  of  the  instrument 
(\  =  1304.200  A,  AA,  =  200  mA).  The  interference  profile  has  been  altered  by  subtracting 
the  average  intensity  of  the  entire  fringe  pattern  from  each  point  in  the  profile.  This  elimi¬ 
nates  much  of  the  DC  component  of  the  fourier  transfonn  as  can  be  seen  in  Figure  21 
(indicated  by  “DC  removed”  in  the  figure  caption).  FFT  spectra  will  now  be  shown  as 
uncalibrated  wavelength  spectra  as  described  above. 


Figure  19  Ideal  Point  Source  Intcrferogram  for  1304.200  A 
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Relotive  Intensity 


Figure  20  Ideal  Point  Source  Interference  Profile  for  1 304.200  A  (DC  removed) 


Figure  21  Ideal  Point  Source  Wavelength  Spectrum  for  1304.200  A  (DC  removed) 
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The  next  series  of  results  (shown  in  Figures  22  and  23)  gives  an  interference 
profile  and  wavelength  spectrum  for  an  incident  wavelength  near  the  nuiximum  allowed 
wavelength  of  the  instrument.  The  FFT  function  produces  an  ambiguous  result  when  the 
wavelength  increases  beyond  the  maximum  allowed  wavelength.  This  ambiguity  is  caused 
by  contamination  of  the  wavelength  peak  due  to  FFT  aliasing  and  is  well  understood.  By 
noting  when  the  aliasing  begins,  the  maximum  wavelength  of  the  instrument  can  be  deter¬ 
mined. 

At  this  point,  sufficient  runs  had  been  conducted  at  known  wavelengths  to 
roughly  calibrate  the  horizontal  axis  of  the  wavelength  spectrum.  The  wavelength  axis 
was  calibrated  by  correlating  the  pixel  position  of  peaks  obtained  on  the  wavelength  spec¬ 
trum  with  the  known  wavelengths  that  produced  them.  The  pixel  corresponding  to  zero 
spatial  frequency  (the  zero  position  on  the  FFT  spectrum)  was  labeled  as  A  linear  rela¬ 
tionship  between  FFT  pixel  spacing  and  wavelength  was  assumed  for  this  project.  A  more 
accurate  calculation  of  the  pixel  spacing  as  a  function  of  wavelength  will  be  required  for 
the  working  prototype  and  flight-qualified  versions  of  the  AMI.  Figure  24  shows  the 
wavelength  spectrum  for  an  ideal  point  source  at  X  -  1304.380  A  with  the  wavelength  axis 
calibration  in  place. 
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Relative  intensity 


Figure  22  Ideal  Point  Source  Interference  Profile  for  1304.380  A 


1304.00  1304.10  1304.20  1304.30  1304.40  1304.50 

Wavelength  (Angstroms) 

Figure  24  Ideal  Point  Source  (Calibrated)  Wavelength  Spectrum  for  1304.380  A 

Several  interference  runs  were  conducted  as  described  above  to  more  accurately 
pinpoint  the  maximum  and  minimum  allowed  values  of  AX.  These  values  were  found  to  be 

AXj^^X  =  400mA  (35) 

and 

AXnjjn  =  3mA.  (36) 

2.  Ideal  Point  Source  with  a  Multiple  Wavelength  Light  Source 

After  the  minimum  and  maximum  allowed  values  of  AX  were  found  and  an 
accurate  wavelength  calibration  scale  introduced,  a  light  source  consisting  of  two  distinct 
wavelengths  was  modeled  in  order  to  determine  the  resolution  of  the  instrument.  The  reso¬ 
lution  here  is  considered  to  be  the  minimum  separation  between  two  wavelengths  such 
that  each  individual  peak  in  the  fourier  domain  is  recognizable.  The  interferogram  result- 
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ing  from  tw'o  equal-intensity  wavelengths  is  shown  in  Figure  25.  Figures  26  and  27  show 
the  resulting  interference  profile  and  wavelength  spectrum.  In  the  Interferogram  and  inter¬ 
ference  profile,  two  separate  spatial  frequencies  are  readily  discernible.  In  the  wavelength 
spectrum  in  Figure  27,  both  wavelength  peaks  are  distinct. 

Figures  28  through  3 1  show  the  wavelength  spectra  resulting  from  pairs  of 
wavelengths  near  130?  800  A  separated  by  gradually  decreasing  wavelength  differences. 
The  wavelengths  became  indistinguishable  when  the  resolution  of  the  instrument  was 
reached.  As  can  be  seen  in  Figure  30,  the  approximate  resolution  is  4.0  mA.  By  conduct¬ 
ing  several  series  of  runs  at  different  wavelengths,  it  was  found  that  the  resolution  is 
roughly  constant  throughout  the  bandpass  of  the  instrument. 


Figure  25  Ideal  Point  Source  Interferogram  for  1303.980  A  and  1303.900  A 
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Relative  intensity 


Figure  26  Ideal  Point  Source  Interference  Profile  for  1303.980  A  and  1303.900  A 
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Figure  27  Ideal  Point  Source  Wavelength  Spectrum  for  1303.980  A  and  1303.900  A 
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Relative  intensity 


Figure  28  Ideal  Point  Source  for  1303.810  and  1303.800  A  (AX  =  10  mA) 


Figure  29  Ideal  Point  Source  for  1303.805  and  1303.800  A  (AX  =  5  mA) 
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Figure  30  Ideal  Point  Source  for  1303.804  and  1303.800  A  (AA,  =  4  inA) 


Wovelength  (Angsl''oms) 


Figure  31  Ideal  Point  Source  for  1303.8035  and  1303.8000  A  (AA,  =  3.5  inA) 
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3.  Offset  Ideal  Point  Source 

Several  interference  runs  were  conducted  using  a  point  source  offset  slightly 
from  the  focus  of  the  parabola.  Offsets  up  to  250  pm  in  the  negative  X  direction  were 
used.  These  runs  were  made  on  the  AMI  version  containing  larger  than  realistic  elements 
so  that  all  rays  exiting  the  mask  of  the  parabola  would  strike  the  detector.  In  this  way,  the 
effect  of  the  offset  could  be  observed  without  the  loss  of  rays  that  would  occur  with 
smaller  components.  It  was  discovered  that  offsets  greater  than  approximately  100  pm 
caused  noticeable  curvature  of  the  fringes.  An  example  of  the  fringe  pattern  caused  by  an 
ideal  point  source  offset  250  pm  from  the  focus  is  shown  in  Figure  32.  Tliis  curvature 
caused  a  deterioration  of  the  interference  profile  since  the  fringes  were  no  longer  parallel 
with  the  vertical  summing  bins  used  by  DART,  making  the  FFT  peak  much  less  pro¬ 
nounced.  The  curvature  in  the  fringe  pattern  will  later  be  used  to  help  explain  the  effect  of 
a  pinhole  on  the  interference  patterns. 


Figure  32  Ideal  Point  Source  (Offset  250  pm)  Interferograin  for  1304.015  A 
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4.  Pinhole  Source 


The  effect  of  a  pinhole  source  on  interference  is  the  eventual  washing  out  of 
fringe  pattents  caused  by  overlapping  fringes  from  multiple  rays  separated  at  the  source. 
Like  the  offset  point  source,  the  pinhole  interference  nms  were  conducted  on  the  large  ele¬ 
ment  AMI  simulation.  Circular  pinholes  were  used  with  diameters  of  100, 250  and  500  pm 
(corresponding  to  commercially  available  sizes).  Figures  33  through  35  show  interfero- 
grams  for  these  pinhole  diameters.  The  100  pm  pinhole  produced  an  interferogram  nearly 
identical  to  that  of  the  ideal  point  source.  Ail  fringes  were  distinct  and  showed  high  con¬ 
trast.  The  interferogram  for  the  250  pm  pinhole  had  a  fringe  contrast  that  was  noticeably 
lower  than  the  100  pm  case.  Fringes  became  almost  conpletely  washed  out  and  indistin¬ 
guishable  using  the  500  pm  source.  Figures  36  and  37  show  the  interference  profile  and 
wavelength  spectrum  for  a  pinhole  diameter  of  500  pm.  Although  the  wavelength  peak  is 
visible  in  Figure  37,  it  is  less  pronounced  than  that  obtained  with  the  ideal  point  source  or 
with  smaller  pinholes. 

While  conducting  this  analysis,  it  became  apparent  that  the  distance  between 
optical  elements  in  the  instrument  design  is  a  major  factor  contributing  to  the  washout 
effect.  The  sensitivity  of  the  instrument  to  fringe  pattern  washout  is  believed  to  be  propor¬ 
tional  to  the  overall  path  length  of  each  beam.  In  the  present  case,  where  the  instrument 
modeled  has  an  especially  large  grating-to-parabola  distance  (60  cm),  the  washout  effect 
was  severe.  Additionally,  distances  between  the  grating  and  the  mirrors  were  larger  than 
necessary  in  this  instrument,  compounding  the  effects  of  fringe  washout.  Although  not 
verified,  it  is  expected  that  shorter  distances  between  the  elements  would  lessen  the  detri¬ 
mental  effects  of  the  size  of  the  pinhole.  Minimizing  the  distance  between  optical  ele¬ 
ments  is  not  expected  to  adversely  affect  instrument  resolution  or  bandpass. 
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It  should  be  pointed  out  that  anomalous  results  were  obtained  when  the  AMI 
was  modeled  with  smaller  more  realistic- sized  elements  (25  mm  diameter  mirrors  with  a 
25x10  mm  detector).  The  fringes  for  this  case  were  strangely  indiscernible,  even  with  the 
100  pm  source.  Although  the  cause  of  this  is  not  fully  understood,  it  is  suspected  that  an 
inadvertent  error  in  the  construct  of  the  parameter  files  contributed  to  the  poor  results 
obtained. 


Figure  33  100  pm  Diameter  Pinhole,  interferogram  for  1304.015  A 
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Figure  35  500  ^im  Diameter  Pinhole,  Interferogram  for  1 304.0 1 5  A 


5.  Slit  Source 


The  final  source  presented  is  the  “real  line  source”  or  slit  source.  This  was 
examined  to  determine  the  effect  on  interference  patterns  from  the  use  of  commercially 
available  slits  at  the  instrument  aperture.  It  was  discovered  that  the  horizontal  dimension 
of  the  slit  produces  instrument  behavior  similar  to  that  of  the  pinhole.  The  vertical  dimen¬ 
sion  of  the  slit  had  little  effect  on  perfomiance.  Therefore,  the  AMI  designer  should  make 
the  choice  between  the  use  of  a  slit  or  pinhole  based  on  intensity  requirements  at  the  detec¬ 
tor.  Slit  or  pinhole  dimensions  should  be  made  based  on  the  trade-off  between  decreased 
fringe  contrast  and  increased  intensity.  Provided  sufficient  light  is  available  at  the  detector, 
maximum  sharpness  of  fringes  can  be  obtained  through  the  use  of  pinholes  with  diameters 
of  100  pm  or  smaller. 

6.  Other  Results 

Other  results  obtained  throughout  the  course  of  this  research  include  the  effect 
of  grating  angle  of  incidence  on  instrument  resolution.  Several  interference  runs  were  per¬ 
formed  using  the  interferometer  portion  of  the  spectrograph  alone  using  a  collimated 
beam.  It  was  discovered  that  as  long  as  the  instrument  parameters  were  accurately  calcu¬ 
lated  for  the  desired  tuned  wavelength,  resolution  was  not  significantly  affected  by  the 
angle  of  incidence.  Therefore,  the  angle  of  incidence  should  be  optimized  based  on  instru¬ 
ment  form,  fit,  and  function  requirements  rather  than  instrument  resolution  desires. 

E.  OXYGEN  1304  A  LINE  PROFILE  INTERFERENCE 
1.  Profile  Modeling 

The  solar  01 1304  A  triplet  profile  was  simulated  using  the  method  of  Gladstone 
(1992)  discussed  above.  The  profile  for  the  1304  A  line  was  modelled  using  the  IDL  pro¬ 
cedure  PROF.PRO  and  the  acconq^anying  function  MODEL  1. PRO,  both  included  in 
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Appendix  A.  As  this  thesis  is  only  concerned  with  a  proof-of-concept,  the  simulation  used 
only  a  representative  set  of  values  for  the  line  integrated  flux,  x-offset,  and  x-dispersion 
from  those  given  by  Gladstone  for  =  1304  A.  The  IDL  procedure  first  converts  x-offset 
and  x-dispersion  from  Doppler  velocity  to  wavelength  using  Equations  (8)  and  (9)  and 
then  calculates  the  flux  as  a  function  of  wavelength  using  Equation  (7).  The  procedure 
then  plots  the  profile  intensity  versus  wavelength.  Additional  IDL  procedures  were  written 
to  write  the  profile  data  to  both  postscript  files  for  later  printing  and  data  files  for  transla¬ 
tion  into  DART.  These  IDL  procedures  are  also  contained  in  Appendix  A. 

A  representative  1304.86  A  line  profile  simulation  is  shown  in  Figure  38  below 
for  tiF  =  3830  ergs  cm'^  sec'  * ,  x^jj  =  10.6  km  sec'  * ,  and  x^jf  =  9..S  km  sec'  * . 


Figure  38  Solar  Oxygen  1304  A  Triplet  Line  Profile  Simulation  (after  Gladstone  (1992)) 
2.  ’  '94  A  Line  Profile  Interference  Results 


The  final  set  of  interference  runs  was  conducted  using  an  incident  light  source 
consisting  of  multiple  wavelengths  to  simulate  the  effects  of  the  atomic  oxygen  triplet  pro 


file  at  1304  A.  As  an  initial  test,  one  side  of  the  emission  line  profile  of  Figure  38  above 
was  modeled  by  creating  an  array  of  wavelength  versus  intensity  values  using  from  five  to 
150  array  elements.  These  data  points  were  entered  into  the  DART  program  as  multiple 
wavelengths  of  varying  intensity.  As  all  relative  intensities  in  DART  must  be  entered  as 
integers,  the  intensity  values  of  the  line  profile  were  scaled  and  truncated.  When  the  wave¬ 
length  interval  between  data  points  was  larger  than  the  instrument  resolution,  the  wave¬ 
length  spectra  from  the  resulting  interference  showed  distinct  peaks  for  each  wavelength. 
Although  there  were  gr^}s  between  the  peaks,  the  outline  of  the  peaks  resembled  the  shs^ 
of  the  line  profile  simulated.  As  more  points  were  added,  nearing  the  instrument  resolu¬ 
tion,  the  size  of  the  g2q)s  decreased  but  there  was  substantial  noise  and  a  large  E>C  compo¬ 
nent  present.  When  the  number  of  points  in  the  model  was  increased  so  that  the  interval 
between  data  points  was  less  than  the  instrument  resolution,  the  shape  of  the  line  profile 
became  indistinct.  This  was  due  to  the  presence  of  a  low  frequency  noise  component.  As 
greater  numbers  of  points  in  the  simulation  were  added,  the  FFT  spectra  did  not  improve. 

In  order  to  minimize  the  effect  of  the  low  frequency  contamination  in  the  fourier 
domain,  the  line  profile  was  offset  firom  Xq.  The  value  selected  for  the  center  of  the  new 
line  profile  was  1303.8  A.  This  offset  had  the  effect  of  modeling  a  new  instrument  with  a 
tuned  wavelength  over  200  mA  higher  than  the  wavelength  of  interest.  With  this  change, 
wavelength  spectra  more  accurately  represented  the  incident  emission  line  profile.  The 
line  profile  was  reproduced  particularly  well  when  the  input  profile  was  compressed  along 
the  wavelength  axis  to  make  it  narrower.  In  this  way,  an  entire  line  profile  could  be  entered 
with  wavelength  spacing  smaller  than  the  pixel  spacing  of  the  instrument,  while  requiring 
that  fewer  than  100  wavelength  values  be  entered  in  DART.  Fewer  wavelengths  at  the 
DART  aperture  tended  to  make  the  application  run  more  quickly  and  with  fewer  errors.  An 
example  of  the  compressed  (and  offset)  line  profile  is  shown  in  Figure  39.  The  resulting 
interference  profile  and  wavelength  spectrum  are  shown  in  Figures  40  and  41,  respec- 


55 


lively.  As  can  be  seen  in  Figure  41,  the  profile  in  the  wavelength  spectrum  closely  resem¬ 
bles  that  of  the  incident  profile.  Of  note  is  the  accuracy  of  the  relative  spacing  of  the  peaks 
and  the  depth  of  the  self-reversal  in  the  reproduced  profile.  The  values  used  in  modeling 
the  profile  displayed  in  Figure  39  are  given  in  Table  III.  For  this  simulation  only  50  wave¬ 
length  bins  were  used.  It  should  be  pointed  out  that  the  quality  of  the  line  profile  recovered 
by  this  technique  (Figure  41)  improves  proportionately  with  the  number  of  wavelength 
bins.  Due  to  limitations  in  the  current  computer  setup,  however,  attempts  to  increase  this 
number  were  not  successful. 


Figure  39  Offset  and  Compressed  Line  Profile  Centered  at  1303.8  A 


Relative  inteosity 


Detector  X  Position  (cm) 


Figure  40  Interference  Profile  for  the  Line  Profile  of  Figure  39 


Figure  41  Wavelength  Spectrum  of  the  Line  Profile  of  Figure  39 
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TABLE  m.  OFFSET  LINE  PROFILE  DATA  POINTS 


Relative 

Intensity 

Wavelength 

(A) 

Relative 

Intensity 

Wavelength 

(A) 

Relative 

Intensity 

Wavelength 

(A) 

1 

1303.7885 

105 

1303.7%5 

95 

1303.8046 

2 

1303.7891 

106 

1303.7970 

87 

1303.8051 

3 

1303.7896 

105 

1303.7975 

77 

1303.8055 

5 

1303.7900 

102 

1303.7980 

65 

1303.8060 

8 

1303.7905 

98 

1303.7985 

54 

1303.8065 

12 

1303.7910 

94 

1303.7991 

43 

1303.8070 

17 

1303.7915 

92 

1303.7994 

33 

1303.8075 

24 

1303.7921 

91 

1303.8000 

25 

1303.8080 

33 

1303.7925 

91 

1303.8004 

18 

1303.8086 

43 

1303.7931 

94 

1303.8010 

12 

1303.8090 

54 

1303.7935 

98 

1303.8016 

8 

1303.80% 

65 

1303.7941 

102 

1303.8020 

5 

1303.8099 

76 

1303.7944 

105 

1303.8026 

3 

1303.8105 

87 

1303.7950 

106 

1303.8030 

2 

1303.8109 

95 

1303.7955 

105 

1303.8036 

1 

1303.8115 

102 

1303.7960 

102 

1303.8040 

0 

1303.8120 
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F.  RESULTS  SUMMARY 


Table  IV  sununarizes  the  numerical  results  obtained  in  conducting  this  research  and 
represents  the  expected  performance  parameters  of  the  AMI.  The  resolving  power  given  is 
calculated  from 

*  =  A  = 

TABLE  IV.  INSTRUMENT  PARAMETERS  BEFORE  CORRECTION 


Parameter 

Result 

max  AX 

400  mA 

min  AX 

3  mA 

Bandwidth 

397  mA 

Resolution 

4mA 

Resolving 

Power 

3x10^ 

Pixel  spacing  is  the  wavelength  spacing  between  elements  in  the  FFT  array.  Alterna¬ 
tively,  it  can  be  thought  of  as  the  smallest  that  will  produce  a  peak  in  the  first  FFT  array 
position  (nearest  to  Xfl).  A  formula  for  pixel  spacing,  given  by  Wallace  (1992)  is 


r' 


AX, 


MIN 


2L 


1  ^0 


(38) 


where  L  is  the  length  of  the  detector.  Substituting  the  values  used  for  this  instrument  into 
Equation  (38)  gives 

AXmin  =  4.0003xl0“^^m^ixel  =4mA/pixel  .  (39) 
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Alternatively,  the  pixel  spacing  can  be  found  experimentally  (or  in  this  case  by  simulatimi) 
by  dividing  the  demonstrated  bandpass  of  the  instrument  by  the  number  of  elements  in  the 
FFT  array.  In  the  case  of  the  AMI,  the  pixel  spacing  is 

The  apparent  discrepancy  between  the  expected  pixel  spacing  and  that  determined  by 
modeling  the  AMI  can  be  explained  by  a  factor-of-2  error  that  was  recently  discovered  in 
the  interference  algorithm  of  the  DART  program  (D.  D.  Qeary,  private  communication, 
1993).  Equation  (23)  gives  the  irradiance  resulting  from  interference  at  a  field  point,  P. 
Factoring  out  the  21  teim,  the  irradiance  is 

Ip  =  21(1  +cos8)  .  (41) 

Through  the  use  of  a  common  trigonometric  identity,  it  becomes 

Ip  =  2Icos^2-  (^2) 

In  its  interference  algorithm,  DART  was  found  to  have  used  the  equation 

Ip  =  2Icos^5.  (43) 

This  introduces  a  factor-of-two  error  in  the  spatial  frequency  of  the  interference  calcula¬ 
tions.  This  error  had  the  effect  of  halving  die  bandwidth  and  the  pixel  resolution  obtained 
for  this  instrument.  Were  this  factor-of-two  error  corrected  for  and  the  interference  runs 
repeated,  the  revised  performance  parameters  would  have  been  those  given  in  Table  V. 
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TABLE  V.  CORRECTED  INSTRUMENT  PARAMETERS 


Parameter 

Result 

max  AX 

800  mA 

min  AX 

3  mA 

Bandwidth 

797  mA 

Resolution 

8  mA 

Resolving 

Power 

1.5  X  10^ 

It  should  be  emphasized  that  the  values  presented  above  in  Tables  IV  and  V  are  based 
on  a  detector  length  of  25  mm  and  a  detector  pixel  spacing  of  60  pm.  The  latter  is  con¬ 
strained  primarily  by  a  limitation  in  the  DART  program.  The  use  of  a  conunercially  avail¬ 
able  one-dimensional  array  detector  with  a  detector  pixel  spacing  of  25  pm  could 
potentially  improve  the  bandwidth  and  resolution. 
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VI.  CONCLUSIONS  AND  RECOMMENDATIONS 


A.  SUMMARY  OF  FINDINGS  AND  CONCLUSIONS 

The  design  of  a  high  resolution  All-reflection  Michelson  Interferometer  for  far  and 
extreme  ultraviolet  atmospheric  spectroscopy  was  studied  using  software  ray  tracing  tech¬ 
niques.  Results  for  various  apertures  and  wavelengths  of  incident  light  were  presented  in 
several  forms.  Interferograms  were  presented  to  show  fringe  patterns  created  at  the  two 
dimensional  detector.  Interference  profiles  displayed  a  plot  of  the  intensity  of  interfering 
light  versus  the  detector  horizontal  axis  position.  Finally,  so-called  wavelength  spectra 
were  produced  by  taking  the  fast  fourier  transform  of  the  interference  profiles. 

An  ideal  point  source  emanating  monochromatic  light  at  the  iqjerture,  produced  high- 
contrast  interference  fringes  on  the  interferogram  and  sharp  peaks  on  the  wavelength  spec¬ 
trum.  Light  of  varied  wavelength  was  used  to  determine  the  minimum  and  maximun> 
allowed  AA.  of  the  instrument  and  a  calibrated  wavelength  scale  was  introduced  for  subse¬ 
quent  wavelength  spectra.  Distinct  peaks  were  observed  on  the  wavelength  spectra  for 
incident  light  consisting  of  two  wavelengths.  Spacing  of  the  incident  wavelengths  was 
varied  to  determine  the  instrument  resolution.  The  expected  bandwidth,  maximum  and 
minimum  allowed  AA,,  resolving  power  and  resolution  of  the  instrument  were  thereby 
determined  and  then  corrected  for  a  factor-of-two  error  discovered  in  the  ray  tracing  soft¬ 
ware.  These  results  were  presented  in  Table  V. 

It  was  found  that  using  a  pinhole  or  slit  as  the  eq)erture  at  the  focus  of  the  off-axis 
paraboloid  tended  to  cause  a  marked  decrease  in  fringe  contrast  compared  with  that  of  an 
ideal  point  source.  It  was  determined  that,  for  the  instrument  modeled  here,  the  maximum 
pinhole  diameter  for  satisfactory  fringe  contrast  was  100  pm.  It  was  also  shown  that  the 
effect  of  offsetting  the  source  of  incident  light  from  the  focus  of  the  collimating  paraboloid 
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is  to  cause  a  curvature  of  the  fringe  pattern.  This  curvature  is  believed  to  be  the  cause  of 
the  decrease  in  fringe  contrast  for  the  wavelength  spectrum  of  the  pinhole  source.  It  is  also 
believed  that  the  sensitivity  of  fringe  contrast  to  pinhole  diameter  could  be  minimized  by 
reducing  the  overall  path  lengths  of  the  instrument.  This  could  be  done,  without  any 
adverse  effect  on  instrument  resolution,  by  reducing  the  spacing  between  optical  elements. 

The  profile  of  the  01  1304  A  emission  was  modeled  on  the  personal  computer  using 
the  method  of  Gladstone  (1992).  The  resulting  profile  was  simulated  on  the  DART  pro¬ 
gram  and  interference  runs  were  conducted.  Fourier  transform  techniques  were  used  to 
determine  the  feasibility  of  retrieving  meaningful  data  from  the  instrument.  It  was  found 
that  the  tuned  frequency  of  the  AMI  should  be  offset  approximately  half  the  bandwidth  of 
the  instrument  from  the  wavelength  of  interest.  This  places  the  profile  in  the  center  of  the 
resulting  wavelength  spectrum.  This  was  effectively  demonstrated  by  modeling  a  scaled 
emission  line  profile  offset  from  Xo.This  profile  was  accurately  retrieved  by  taking  the  fou- 
rier  transform  of  the  interference  profile  obtained  from  the  AMI.  Finally,  the  suspected 
factor-of-two  error  in  the  DART  application  was  confirmed  by  examination  of  the  pixel 
spacing  of  the  instrument. 

The  results  presented  herein  have  validated  the  overall  design  concept  of  the  Plane 
Grating  All-reflection  Michelson  Interferometer  as  a  viable  means  of  examining  atomic 
oxygen  emission  line  profiles.  It  is  believed  that  with  further  development,  this  instrument 
has  potential  to  play  an  important  role  in  determining  atomic  oxygen  distributions  in  the 
atmosphere  and  to  become  a  catalyst  in  the  further  understanding  of  the  physics  of  the 
atmosphere. 

B.  RECOMMENDATIONS  FOR  FURTHER  STUDY 

The  research  presented  in  this  thesis  was  conducted  using  purely  simulation  and 
modeling  techniques  on  the  personal  computer.  Further  work  must  now  be  done  to  build 
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and  test  a  working  prototype  of  this  instrument  in  the  visible  and  ultraviolet  regions  of  the 
electromagnetic  spectrum,  continuing  on  the  work  performed  in  the  laboratory  by  Wallace 
(1992)  and  Risley  (1993).  Of  particular  interest  is  experimentation  with  shorter  coherence 
length  sources  and  the  testing  of  the  fourier  transform  techniques  used  here  on  actual 
exfyerimental  data. 

Further  work  in  the  simulation  area  includes  revising  the  DART  application  to  ensure 
that  a  sigitiiicantly  large  number  of  wavelengths  may  be  entered  at  the  DART  aperture.  A 
revision  to  allow  the  user  to  select  the  number  of  bins  under  which  the  interference  data  is 
summed  is  also  required.  By  so  doing,  equivalent  resolutions  can  be  achieved  by  simula¬ 
tion  and  by  experimentation  usmp  a  1024  pixel  detector.  Additional  work  in  the  immediate 
area  of  this  thesis  can  include:  ( 1 )  devising  a  method  for  calculating  and  then  removing  the 
increase  in  intensity  of  the  data  with  x-axis  position  caused  by  the  off-axis  parabola;  (2) 
designing  an  AMI  that  is  tuned  to  a  wavelength  offset  from  the  desired  so  that  wave¬ 
lengths  on  either  side  of  Xo  can  be  distinguished,  providing  a  two-sided  representation  of 
the  line  profile  of  interest;  (3)  investigating  the  source  and  removal  of  low  frequency  con¬ 
tamination  in  wavelength  spectra;  (4)  continuing  the  examination  of  the  effects  on  instru¬ 
ment  resolution  from  pinhole  sources  and  slits  at  the  aperture;  (5)  redesigning  the 
instrument  to  reduce  the  spacing  between  optical  elements  and  therefore  minimize  the  det¬ 
rimental  effects  of  pinhole  size  on  fringe  contrast;  and  (6)  designing^the  entrance  optics  to 
narrow  the  bandwidth  of  the  light  incident  on  the  AMI ’s  grating  so  that  the  numerous 
orders  caused  by  the  continuum  of  wavelengths  at  the  instrument’s  entrance  slit  will  not 
unintentionally  strike  the  detector  and  corrupt  the  data  obtained. 
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APPENDIX  A 


IDL  FUNCTION  AND  PROCEDURE  LISTINGS 


IDL  Function  MODELl.PRO 

This  function  provides  the  equation  used  to  model  the  1304  A  line  profile. 

FUNCTION  MODEL1.FLUXLDISXOFF.LO.L 
02.998E8 

LEFr*FLUX/(2*SQRT(  !PI)*((C*LD1S/L0)-C)) 

RIGHT=EXP(-((L-LOFFy(LDIS-LO)y'2)+EXP(-((L-(2*LO)+LOFFy(LDIS-LO)y'2) 

RETURN.(LEFT*RIGHT) 

END 


IDL  Procedure  PROF.PRO 


;  A  program  to  plot  the  oxygen  1304  A  line  profile  to  screen 


@MODELl 

FLUX*3830.0 

XDIS»10.6 

XOFF»=9.5 

LAMBDACTR=1304.0 

NUMPOINTS-40 

WIDTH*2.0E-10 

C=2.998E8 

LAMBDA0=LAMBDACrR*  1.E- 10 

XDIS=XDIS*1000 

XOFF^XOFF*1000 


;IN  KM/SEC 
;IN  KM/SEC 
;IN  ANGSTROMS 

;RES0LUT10N  IS  WIDTH/NUMPOINTS 

;BANDWIDTH  IN  METERS 

;INM/S 

;LAMBDA  IN  METERS 
:GIVES  X’s  IN  METERS/SEC 


LAMBDADIS=(XDIS+C)*LAMBDAO/C 

LAMBDAOFF=(XOFF+C)*LAMBDAO/C 

HALFWIDTH=WIDTH/2.0 

LAMBDA=FLTARR(NUMPOINTS) 

PROFILE»FLTARR(NUMPOINTS) 

LAMBDA=((FINDGEN(NUMPOINTS)*WIDTH/NUMPOINTS)+LAMBDAO-HALFWIDTH) 
PROnLE=MODELl(FLUXiAMBDADISLAMBDAOFF,LAMBDAO, LAMBDA) 


PLOT,LAMBDA*1.E10JPROFILE,XRANGE=(1303.5.I304.5]3ACK=255,COLOR«4),$ 
TTTLEs’SOLAR  oxygen  1304  ANGSTROM  TRIPLET  LINE  PROHLE’.S 
XTTTLEs’WAVELENGTH  (ANGSTR0MS)’.$ 

YnTLE=’COUNT  RATE  (COUNTS/S/MILLIANGSTROM)’ 

END 
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IDL  Procedure  PROFFIL.PRO 


;A  prognun  to  wiite  the  oxygen  1304  A  profile  ditta  to  a  file  for  inclusion  later  in  a  DART  parameter  file 


@M0DEL1 

FLUX=3830.0 

XDIS*10.6 

XOFP=9.5 

LAMBDACrR*1304.0 

NUMPOINTS*40 

WIDTH»0.3E-10 

C-2.998E8 

LAMBDAO-LAMBDACTR*  1  E- 10 

XDIS»XDIS*1000 

XOFF»XOFF*1000 


;IN  KM/SEC 
;IN  KM/SEC 
;IN  ANGSTROMS 

;RESOLUTION  IS  WWTH/NVMPOINTS 
;BANDWIDTH  IN  METERS 
:1NM/S 

■.LAMBDA  IN  METERS 
;GIVES  X’s  IN  METERS/SEC 


LAMBDADIS»(XDIS+C)*LAMBDAO/C 

LAMBDAOFF«(XOFF+C)*LAMBDAO/C 

HALFWIDTH*WIDTH/2.0 

LAMBDA«FLTARR(NUMPOINTS) 

PROFILE*FLTARR(NUMPOINTS) 

LAMBDA-((FINDGEN(NUMPOINTS)*WIDTH/NIJMPOINTS)+LAMBDAO-HALFWIDTH) 
PROFILE»MODELl(FLUXaAMBDADISAAMBDAOFF.LAMBDAO.LAMBDA) 
PROFILE2»FIX(PROFILE*  1000.) 

OPENW.l  ,’OXYGEN.DAT’ 

FOR  J-0.NUMPOINTS.1  DO  PRINTF.1JFORMAT*=’{IJP10.4)’,PROFILE2(J).LAMBDA(J)*1.E10 

CLOSE,  1 

END 
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IDL  Procedure  PROFPRT.PRO 

;Ptooedure  to  plot  the  Oxygen  1304  A  line  profile  to  a  postsciipt  file  for  printing  as  a  figure 


@MODELl 

FLUX-3830.0 

XDIS-10.6 

XOFF-9.5 

LAMBDACTR-1304.0 

NUMPOINTS-39 

WIDTH-2.0E-10 

C-2.998E8 


;IN  KM/SEC 
:IN  KM/SEC 
.IN  ANGSTROMS 

;RESOLUnON  IS  WIDTH/NUMPOINTS 
;BANDWIDTH  IN  METERS 
ON  M/S 


LAMBDAO-LAMBDACTR^l.E'lO  O^AMBDA  IN  METERS 

XDIS-XD1S*1000  ;GIVES  X*s  IN  METERS/SEC 

XOFF»XOFF*1000 

LAMBDADIS=(XDIS+C)*LAMBDAO/C 

LAMBDAOFF»(XOFF+C)*LAMBDAO/C 

HALFWIDTH-WIDTH/2.0 

LAMBDA-FLTARR(NUMPOlNTS) 

PROFILE«FLTARR(NUMPOINTS) 

LAMBDA-((FINDGEN(NUMPOINTS)*WIDTH/NUMPOINTS)+LAMBDAO-HALFWIDTH) 

PROFILE-MODEL1(FLUXXAMBDADISAAMBDAOFF,LAMBDAO.LAMBDA) 


SET.PLOT.’PS’ 

DEVICE,FILENAME»’HGURE.PS’yiNCHESyP0RTRAIT,SCALE-.72 
PLOT.LAMBDA*l.E10J'ROFILE.XRANGE-tl303.5.1304.5]3ACK=255,COLOR*0.$ 
TITLE-’SOLAR  OXYGEN  1304  ANGSTROM  TRIPLET  LINE  PROFILE’, $ 
XnTLE-’WAVELENGTH  (ANGSTROMS)’,$ 

YTITLE-’COUNT  RATE  (COUNTS/S/MILLIANGSIROM)’ 

DEVICE./CLOSE 

SET.PLOT.’WIN’ 

END 
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IDL  Procedure  INPLOT.PRO 


;This  is  an  IDL  procedure  that  will  read  a  file  produced  by  DART 
consisting  of 400  vertical  bins  of  interfereiK:e  intensities 
along  the  x  direction,  plot  the  resulting  data,  take  a  fourier 
transform  of  the  figure,  then  plot  the  results  all  to  the  screen. 

XAX»FLTARR(400)  :PLOT  X-AXIS  DATA  ARRAY 

XPOS«FLTARR(400)  ;POSrnON  ALONG  DETECTOR  X-AXIS 

INTEN=FLTARR(400)  :INTENSITY  (Y-AXIS)  ARRAY 

PEAKSsO 
JUNK»STRING(0) 

FILENAMEsSTRINGfO) 

PRINT, 'ENTER  THE  NAME  OF  THE  INTERFERENCE  FILE  ‘ 

READJTLENAME 
OPENR,l, FILENAME 

READF,14UNK  ;USE  UP 'DART  INTERFERENCE  FILE’ TEXT 

FOR  1>0,399  DO  BEGIN 
READF,IXY 
XPOS(I)=X 
INTEN(I)»Y 
ENDFOR 
CLOSE,  1 

XMAX«MAX(XPOS) 

XMIN>MIN(XPOS) 

:SUBTRACT  THE  AVERAGE  VALUE  FROM  THE  INTENSITY 
lARRAYTO  MAKE  IT  SYMMETRIC  ABOUT  THE  X  AXIS. 
INTEN«INTEN-(TOTAL(INTEN)/400.) 

;USE  ‘REVERSE’  IF  PLOTTING  LEFT  SIDE  OF  LAMBDA  0 
;USE  REGULAR  XAX  IF  PLOTTING  RIGHT  SIDE  OF  LAMBDA  0 
:XAX=REVERSE(1304.-{.00207*FINDGEN(200))) 

XAX=1304.+(.00207*FINDGEN(200)) 

PLOT,XPOS,INTEN,XRANGE*[XMIN,XMAX]3ACK=255,COLOR=l,$ 
TTTLE=’INTERFERENCE  PATTERN’,XnTLE*’Position  on  Detector  (cm)’,$ 

YTITLE® ’Relative  Intensity’ 

;THE  FOLLOWING  SECHON  NOT  USED  IN  FINAL  PRODUCT  BUT  ALLOWS 
;USER  TO  EXAMINE  THE  FIRST  PLOT  BEFORE  CALLING  UP  THE  SECOND 
PRINT, ’ENTER  THE  NUMBER  OF  FRINGES  (PEAKS)’ 

READ3EAKS 

;L=K(XMAX-XMIN)/PEAKS)/100. 

;DELLAMBDA*(1304.E-10/(2.*L))*SQRT((l./(6.E6)y'2-(13(>4.E-10/2.y'2) 
;PRINTJ>ELLAMBDA*1.E10,’ ANGSTROMS’ 

FORT«ABS(FFT(INTEN,-l))  -.FAST  FOURIER  TRANSFORM 

PLOT,XAX,FORT(0;199)3ACK=255,COLOR=l,TrTLE=’FREQUENCYSPECTRUM’.$ 
XTTTLE=’Wavelength  (Angstroms)’,YliTLE=’Relative  Intensity’, $ 
XRANGE=[1304.0,1304.4] 

END 
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IDL  Procedure  INPRT.PRO 


This  procedure  plots  the  interference  profile  and  the  FFT  to  a  postcript  file  but  doesn’t 
subtract  off  the  average  intensity  to  remove  the  unwanted  DC  component  nor  does  it  scale 
the  wavelength  axis.  The  FFT  plots  produced  here  are  double  sided. 

;Tbis  is  an  IDL  program  that  will  read  a  file  produced  by  DART 
consisting  of 400  bins  of  interference  intensities  summed 
along  the  y  direction,  plot  the  resulting  data,  tidre  a  fourier 
transform,  then  plot  die  results  all  to  a  postscript  file. 

XPOS«FLTARR(400) 

INTEN»FLTARR(400) 

JUNK«STRING(0) 

FILENAME-STRING(0) 

PRINT,’ENTER  THE  NAME  OF  THE  INTERFERENCE  FILE  ‘ 

READ,FILENAME 
OPENR,l, FILENAME 

READRIJUNK  ;USE  UP ‘DART  INTERFERENCE  FILE’ TEXT 

FOR  I»0,399  DO  BEGIN 
READF,1JC,Y 
XPOS(I)-X 
INTEN(I)»Y 
ENDFOR 
CLOSE,  1 

XMAX-MAX(XPOS) 

XMIN*MIN(XPOS) 

SEr_PLOT,’PS’ 

DEVlCE,HLENAME»’HG2.PS’7INCHES7PORTRAIT,SCALE».72 

PLOT,XPOS,INTEN,XRANGE=t-2.,2.]3ACK»255,COLOR»l,$ 
titles’  INTENSITY  VS.  DETECTOR  X  POSITION*  ,$ 

X 1 1 1  iJEs ’Detector  X  Position  (cm)’,YnTLE=’Relative  Imensity’ 

FORTsABS(FPT(INTEN,- 1 )) 

PLOT,FORT3ACK=255,COLOR*l,TITLE=’  FOURIER  TRANSFORM  PLOT‘,$ 

XTTTLEs’FFT  Solution  Array  Element  Number’,YnTLE= ’Relative  Intensity’,! 
XRANGEslO.00,400] 

PLOT,FORT3ACKs255,COLOR=1  .TITLEs’IDEAL  PS  FFT  PLOT’,! 

XTITLEs’FFT  Solution  Array  Position’,Y  ll  liJE=’Relative  Intensity’,! 

XRANGEs[180,200] 

DEVICE,/CLOSE 

SET.PLOX’WIN’ 

END 
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IDL  Procedure  IPRT.PRO 

This  procedure  plots  the  interference  profile  and  the  FFT  to  a  postcript  file  and  then  sub¬ 
tracts  off  the  average  intensity  to  remove  the  unwanted  DC  component.  It  also  introduces 
the  wavelength  axis  calibration  and  plots  the  wavelength  spectrum. 


;This  is  an  IDL  program  that  will  read  a  file  produced  by  DART 
coosistiag  of  400  bins  of  interference  intensities  summed 
along  the  y  direction,  plot  the  resulting  data,  take  a  fourier 
tmsfoim  ,  then  plot  the  results  all  to  a  postscript  file. 

XAXsPLTARR(400) 

XPOS»FLTARR(400) 

INTEN*FLTARR(400) 

JlJNK=STRING(0) 

FILENAME»STRING<0) 

PRINT, ’ENTER  THE  NAME  OF  THE  INTERFERENCE  FILE  ‘ 

READJ^ILENAME 

OPENR,lJ!ILENAME 

READEIJUNK  :USE  UP ’DART  INTERFERENCE  PILE’ TEXT 

FOR  1-0,399  DO  BEGIN 
READF,1,X,Y 
XPOS(I)*X 
INTENa)*Y 
ENDFOR 
CLOSE,  I 

XMAX«MAX(XPOS) 

XMIN«MIN(XPOS) 

SET_PLOT,’PS’ 

DEVICE,FILENAME=’t7.PS’yiNCHESyPORTRAIT,scale=.72 

INTEN=INTEN-(TOTAL(INTEN)/400.) 

;XAX-1304.+(.00207*FINDGEN(200)) 

XAX=REVERSE(1304..(.00207*FINDGEN(200))) 

PLOT,XPOS,INTEN,XRANGE=[-2.,2.]3ACK-255,COLOR=l,$ 

TITLE-’  INTENSITY  VS.  DETECTOR  X  POSITION’  ,$ 

XTTTLE-’Detector  X  Position  (cm)’,YTITLE=’Relalive  Intensity’ 

FORT=ABS(FFr(INTEN,- 1 )) 

PLOT,XAX,FORT(200:*)3ACK=255,COLOR»1,TITLE=’  FFT  PLOT‘,$ 

Xiiii.Es’ Wavelength  ( Angstroms)’ ,Y  ill  LE- ’Relative  Intensity’,$ 

XRANGEs{1303.6,1304.01 

PLOT.XAX,FORT(200:*)3ACK=235,COLORsl,TrrLE=’  FFT  PLOT’,$ 

X’iiii.Es’ Wavelength  ( Angstroms  )’,YTli'L£s ’Relative  lntensity’,$ 
XRANGEs[1303.75,1303.85] 

DEVICE,/CLOSE 

SET.PLOT’WIN’ 

END 
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APPENDIX  B 


DART  PARAMETER  FILES 


Parameter  File  OORDER.TXT 

This  file  contains  the  parameters  for  a  typical  run  for  the  zero  order  path  of  the  original 
interferometer. 

04/16/1993  18:13:06 

The  analysis  wavelength  is  1304.0001  angstroms. 

Aperture  data: 

Center  of  Opening  (x,y,z):  [  0.00000, 0.00000, 0.00000] 

Aperture  dimensions  (x,y):  [  2.50000, 2.S0000] 

Points  across  opening  (x,y):  ( 10, 10] 

Angular  beam  spread  (x,y):  [  0.000, 0.000] 

Number  of  angles  in  qnead  (x,y);  [1,1] 

Center  angle  of  beam  rays  (x,y):  [  0.000, 0.000] 

No  offset  aiigles  are  traced. 

Aperture  shape  is  rectangular. 

Beam  spread  sh^  is  rectangular. 

There  are  4  elements: 
xy  z 

The  Center  of  Element  1  coordinate  System  is  0.000000  0.000000  40.000000 
Euler  Angles  in  Degrees:  0.000000  147.000000  0.000000 
Plane  Grating  with  60000.000  grooves  per  cm  atxi  order  0. 

Element  mask  is  rectangular  with  width  5.000  and  height  5.000  cm. 

X  y  z 

The  Center  of  Element  2  coordinate  System  is  18.271000  0.000000  31.865000 
Euler  Angles  in  Degrees:  0.000000  -66.000000  0.000000 
Plarre  Mirror 

Element  mask  is  circular  with  width  5.000  and  height  5.000  cm. 

X  y  z 

The  Center  of  Element  3  coordinate  System  is  0.000000  0.000000  40.000000 
Euler  Angles  in  Degrees:  0.000000 147.000000  0.000000 
PlaiK  Grating  with  60000.000  grooves  per  cm  atxl  order  1 . 

Element  mask  is  rectangular  with  width  5.000  and  height  5.000  cm. 
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xyz 

The  Center  of  Element  4  coordinate  System  is  18.21 1000  0.000000  22.872000 
Euler  Angles  in  Degrees:  0.000000  313.246000  0.000000 
Plane  Surface 

Element  mask  is  rectangular  with  width  3.000  and  height  5.000  cm. 

There  were  100  rays  from  all  points  and  angles. 

100  rays  hit  the  detector. 

0  failed  to  hit  first  element. 

0  were  stopped  after  the  first. 


Parameter  File  -10RDER.TXT 

This  file  contains  the  parameters  for  a  typical  run  for  the  minus  first  order  path  of  the  orig¬ 
inal  interferometer. 

03/22/1993  11:04:47 

The  aiudysis  wavelength  is  1305.0000  angstroms. 

Aperture  data: 

Center  of  Opening  (x,y^):  [  0.00000, 0.00000, 0.00000] 

^jerture  diinensions  (x,y):  [  2.50000, 2.50000] 

Points  across  opening  (x,y):  [  100, 100] 

Angular  beam  spread  (x.y):  [  0.000, 0.000] 

Number  of  angles  in  ^read  (x,y):  [  1, 1] 

Center  angle  of  beam  rays  (x,y):  [  0.000, 0.000] 

No  oCEset  angles  are  traced. 

Aperture  shape  is  rectangular. 

Beam  spread  shape  is  rectangular. 

There  are  4  elements: 

X  y  z 

The  Center  of  Element  1  coordinate  System  is  0.000000  0.000000  40.000000 
Euler  t\ngles  in  Degrees:  0.000000  147.000000  0.000000 
Plane  Grating  with  60000.000  grooves  per  cm  aixl  order  -1 . 

Element  mask  is  rectangular  with  width  5.000  and  height  5.000  cm. 

X  y  z 

The  Center  of  Element  2  coordinate  System  is  6.592000  0.000000  21.118000 
Euler  Angles  in  Degrees:  0.000000  340.754000  0.000000 
Plane  Mirror 

Element  mask  is  circular  with  width  5.000  aixl  height  S.OOO  cm. 

X  y  z 

The  Center  of  Element  3  coordinate  System  is  0.000000  0.000000  40.000000 
Euler  Angles  in  Degrees:  0.000000  147.000000  0.000000 
Plane  Grating  with  60000.000  grooves  per  cm  and  order  0. 

Element  mask  is  rectangular  with  width  5.000  and  height  5.000  cm. 


72 


xyz 

The  Center  of  Element  4  coordinate  System  is  1 8.2 1 1000  0.000000  22.872000 
Eukr  Angles  in  Degrees:  0.000000  3 13.246000  0.000000 
Plane  Surface 

Element  mask  is  rectangular  with  width  S.OOO  and  height  S.OOO  cm. 
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Parameter  file  PAR4.TXT 


This  file  contains  the  ofif-axis  parabolic  reflector  element. 

06/03/1993  18:27:20 

The  analysis  wavelength  is  3000.0000  angstroms. 

Aperture  dau: 

Center  of  Opening  (x,y^):  [  0.00000, 0.00000, 0.00000] 

^wituie  dimensions  (x,y):  [  0.01000, 0.01000] 

Points  across  opening  (x.y):  [  1, 1] 

Angular  beam  spread  (x,y):  [  43.000, 43.000] 

Number  of  angles  in  ^nead  (x,y):  [  10, 10] 

Center  angle  of  beam  rays  (x,y):  [  3.000, 0.000] 

No  offset  angles  are  traced. 

Aperture  is  adjusted  to  an  elliptical  shape. 

Beam  spread  shape  is  rectangular. 

There  ate  3  elements: 

X  y  z 

The  Center  of  Element  1  cootdiruie  System  is  -3.970000  0.000000  0.000000 
Eukr  Angles  in  Degrees:  0.000000  90.000000  0.000000 

Concave  Parabrdic  Mirror  with  x-axis  focal  length  «  3.970  and  y-axis  focal  length  k  3.970 
Element  mask  is  circular  with  width  31.000  and  height  31.000  cm. 

xyz 

The  Center  of  Element  2  coordinate  System  is  3.000000  0.000000  11.940000 
Euler  Angles  in  Degrees:  0.000000  90.000000  0.000000 
Plane  Surface 

Element  mask  is  circular  with  width  6.330  and  height  6.330  cm. 
xyz 

The  Center  of  Element  3  coordinate  System  is  20.000000  0.000000  10.000000 
Euler  Angles  in  Degrees:  0.000000  90.000000  0.000000 
Plane  Surface 

Element  mask  is  circular  with  width  20.000  and  height  20.000  cm. 


There  were  100  rays  from  all  points  atxl  angles. 
27  rays  hit  the  detector. 

34  failed  to  hit  first  element 
39  were  stopped  after  the  first. 
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Parameter  file  IPSO.TXT 


This  file  describes  the  zero  order  path  of  the  entire  All-reflection  Micl^lson  Interferometer 
as  simulated  in  this  thesis. 

05/12/1993  09:14:02 

The  analysis  wavelength  is  1303.6300  angstroms. 

Aperture  data: 

Center  of  Opening  (x,y,z):  [  0.00000, 0.00000, 0.00000} 

Aperture  dimensions  (x,y):  [  0.00000, 0.00000] 

Points  across  opening  (x,y):  [  1, 1] 

Angular  beam  spread  (x,y):  [  12.000, 12.000] 

Number  of  angles  in  spread  (x,y):  [  4, 4] 

Ceitter  angle  of  beam  rays  (x,y):  [  -1.000, 0.000] 

No  offset  angles  are  traced. 

Aperture  is  adjusted  to  an  elliptical  shape. 

Beam  spread  is  adjusted  to  an  elliptical  slu^ie. 

There  are  6  elements: 
xy  z 

The  Center  of  Element  1  coordiiutte  System  is  -5.970000  0.000000  0.000000 
Euler  Angles  in  Degrees:  0.000000  90.000000  0.000000 

Concave  Parabolic  Mirror  with  x-axis  focal  length  s  5.970  and  y-axis  focal  length  =  5.970 
Element  mask  is  circular  with  width  3 1.000  and  height  31 .000  cm. 

X  y  z 

The  Center  of  Element  2  coordinate  System  is  5.000000  0.000000  11.940000 
Euler  Angles  in  Degrees:  0.000000  90.000000  0.000000 
Plane  Surface 

Element  mask  is  circular  with  width  6.350  and  height  6.350  cm. 

X  y  z 

The  Center  of  Element  3  coordinate  System  is  60.0000000.000000  11.940000 
Euler  Angles  in  Degrees:  0.000000  237.000000  0.000000 
Plane  Grating  with  60000.000  grooves  per  cm  and  order  0. 

Element  mask  is  rectangular  with  width  5.000  and  heigttt  5.000  cm. 

xy  z 

The  Center  of  Element  4  coordinate  System  is  51.8652700.000000  -6.330910 
Euler  Angles  in  Degrees:  0.000000  24.000000  0.000000 
Plane  Mirror 

Element  mask  is  circular  with  width  5.000  and  height  5.000  cm. 

X  y  z 

The  Center  of  Element  5  coordinate  System  is  60.000000  0.000000  1 1 .940000 
Euler  Angles  in  Degrees:  0.000000  237.000000  0.000000 
Plane  Grating  with  60000.000  grooves  per  cm  and  order  1 . 

Element  mask  is  rectangular  with  width  5.000  and  height  5.000  cm. 
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X  y  z 

The  Center  of  Element  6  coordinate  System  is  42.872000  0.000000  -6.27 1000 
Eukr  Angles  in  Degrees:  0.000000  43.24SS70  0.000000 
Plane  Surface 

Element  mask  is  rectangular  with  width  2.S00  and  height  1 .000  cm. 

There  were  12  rays  from  all  points  and  angles. 

2  rays  hit  the  detector. 

0  failed  to  hit  first  element. 

10  were  stopped  after  the  first. 


Parameter  file  IPS-l.TXT 

This  file  describes  the  minus  first  order  path  of  the  entire  All-reflection  Mkhelson  Interfer¬ 
ometer  as  simulated  in  this  thesis. 

05/26/1993  15:32:58 

Ideal  Point  Source  -  Tuned  to  1304.000  Angstroms  -1  order  ^ath 
The  analysis  wavelength  is  1303.6500  angstroms. 

Aperture  data: 

Center  of  Opening  (x,y,z):  [  0.00000, 0.00000, 0.00000] 

Aperture  dimensions  (x,y):  [  0.00000, 0.00000] 

Points  across  opening  (x,y):  [  1, 1] 

Angular  beam  ^read  (x,y):  [  12.000, 12.000] 

Number  of  angles  in  spread  (x.y):  [  5, 5] 

Center  angle  of  beam  rays  (x,y):  [  -1.000, 0.000] 

No  offset  angles  are  traced. 

Aperture  is  adjusted  to  an  elliptical  shape. 

Beam  spread  is  adjusted  to  an  elliptical  stuqre. 

There  are  6  elements: 

X  y  z 

The  Center  of  Element  1  coordinate  System  is  -5.970000  0.000000  0.000000 
Euler  Angles  in  Degrees:  0.000000  90.000000  0.000000 

Concave  Parabtdic  Mirror  with  x-axis  focal  length  s  5.970  and  y-axis  focal  length  >  5.970 
Element  mask  is  circular  with  width  3 1.000  and  height  31.000  cm. 

X  y  z 

The  Center  of  Element  2  coordinate  System  is  5.000000  0.000000  11.940000 
Euler  Angles  in  Degrees:  0.000000  90.000000  0.000000 
Plane  Surface 

Element  mask  is  circular  with  width  6.350  and  height  6.350  cm. 

X  y  z 

The  Center  of  Element  3  coordinate  System  is  60.000000  0.000000  1 1 .940000 
Euler  Angles  in  Degrees:  0.000000  237.000000  0.000000 
Plane  Grating  with  60000.000  grooves  per  cm  and  order  -1. 

Element  mask  is  rectangular  with  width  5.000  and  height  5.000  cm. 


76 


ny  z 

The  Center  of  Element  4  coordinate  System  is  41.117710  0.000000  5.347650 
Euler  Angles  in  Degrees:  0.000000  70.754430  0.000000 
Plane  Mirror 

Element  mask  is  circular  with  width  5.000  and  height  5.000  cm. 
xy  z 

The  Center  of  Element  5  coordinate  System  is  60.000000  0.000000  1 1 .940000 
Euler  Angles  in  Degrees:  0.000000  237.000000  0.000000 
Plane  Grating  with  60000.000  grooves  per  cm  and  order  0. 

Element  mask  is  rectangular  with  width  5.000  and  hei^t  5.000  cm. 

X  y  z 

The  Center  of  Element  6  coordinate  System  is  42.872000  0.000000  -6.27 1000 
Euler  Angles  in  Degrees:  0.000000  43.245570  0.000000 
Plane  Surface 

Element  mask  is  rectangular  with  width  2.500  and  height  1 .000  cm. 


There  were  23  rays  from  all  points  and  angles. 
6  rays  hit  the  detector. 

0  failed  to  hit  first  element. 

17  were  stopped  after  the  first. 
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